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SUMMARY 

A logic  for  an  NDE  methodology  to  quantitatively  evaluate 
the  integrity  of  structural  elements  is  explored.  The  methodology 
is  based  on  a simplified  analysis  to  predict  flaw  criticality  and 
growth,  an  experimental  program  to  support  and  verify  the  analy- 
sis, and  an  NDE  technique  to  use  NDI  measurements  in  liaison  with 
the  analysis  to  determine  flaw  criticality.  Typical  flaws  are 
identified  and  static  analyses  developed  or  presented  to  predict 
their  growth.  A "wearout"  fatigue  concept  is  incorporated  into 
the  static/fatigue  model  and  the  concept  of  an  instability  boun- 
dary established  for  every  flaw  size  and  location.  An  experimen- 
tal program  is  outlined  to  support  the  analysis  and  to  measure 
the  specific  adhesive  surface  fracture  energy  quantitatively. 
Acoustic  monitoring  of  the  specimens  to  detect  damage  growth  and 
a corresponding  vibration  analysis  are  carried  out  as  a potential 
NDT  and  for  correlation  studies.  Ultrasonic  "C"  scans  are  used 
to  ensure  quality  control  and  to  determine  approximate  damage 
growth  rates.  The  analytical  and  experimental  results  are 
incorporated  into  an  NDE  technique  to  quantitatively  determine 
a measure  of  the  residual  lifetime  and  strength  of  the  structural 
component.  A comparison  of  these  with  a quantified  criterion 
for  damage  tolerance  would  establish  flaw  criticality  through 
the  need,  or  lack  of  it,  for  mandatory  repair. 
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INTRODUCTION 

In  evaluating  the  integrity  of  metallic  structural  ele- 
ments, nondestructive  inspection  has  focused  upon  detection 
of  both  surface  and  subsurface  flaws  and  determination  of  an 
effective  flaw  size.  Given  the  flaw  size,  classical  frac- 
ture mechanics  technology  has  been  employed  to  predict  its 
rate  of  growth  and  the  critical  length  which  corresponds  to 
failure  of  the  element.  Inherent  in  this  approach  is  the 
assumption  that  microscopic  flaws  and  imperfections  grow  into 
macroscopic  cracks  when  the  structural  element  containing 
the  flaw  is  subjected  to  a fatigue  loading  and/or  an  adverse 
environment.  In  addition,  the  macroscopic  isotropy  of  metal- 
lic materials  meant  that  generally  initial  cracks  propagated 
due  to  crack-tip  tensile  stresses  normal  to  the  crack  plane. 
For  contemporary  metallic  materials  and  structures,  tech- 
nology sufficient  to  detect  and  assess  the  criticality  of 
flaws  upon  component  life  ha^s  been  developed.  Of  course, 
there  continues  to  be  a need  to  refine  nondestructive  evalua- 
tion techniques  and  crack  growth  models. 

Fiber-reinforced  composite  materials  are  more  complex 
at  the  macroscopic  level  than  the  contemporary,  metallic 
materials.  These  oaterials  therefore  present  new  problems 
and  challenges,  both  in  flaw  detection  and  in  assessing  the 
influence  of  flaw  characteristics  upon  strength.  Also, 
determination  of  the  influence  of  a given  flaw  geometry  upon 
strength  and  stiffness  is  considerably  more  complex  than  the 
classical  fracture  mechanics  approach  due  to  the  anisotropy 
of  the  fiber-reinforced  composite  materials. 

There  has  been  a recent  interest  in  the  quantitative 
evaluation  of  the  damage  tolerance  of  composite  structures. 
Given  a composite  structural  element,  it  is  desired  to  relate, 
quantitatively,  the  existent  in  situ  flaw  size  and  the  resid- 
ual lifetime  and  the  residual  strength  of  the  structural  ele- 
ment. This  will  enable  definition  of  regular  inspection  in- 
tervals for  the  structural  element  during  its  service  to 
determine  the  criticality  of  the  existent  flaws. 
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Various  techniques  have  been  used  to  assure  the  suit- 
ability of  a structure  for  its  intended  function.  Proof- 
testing to  a load  level  in  excess  of  any  expected  service 
load  is  a possible  approach.  However,  although  proof -testing 
verifies  the  initial  load-bearing  capacity  of  the  structural 
element,  it  cannot  be  used  to  predict  the  growth  of  existent 
flaws  with  fatigue  loading  and  the  effect  of  this  growth  on 
the  residual  strength  and  lifetime  of  the  element.  Spectrum 
fatigue  loading  is  an  alternate  approach  which  would  provide 
the  required  lifetime  data,  but  in  many  cases  is  prohibitive 
cost-wise. 

The  above  drawbacks  in  the  conventional  experimental 
procedure  indicate  the  potential  advantages  of  a new  method- 
ology which  combines  analysis  with  a limited  number  of  experi- 
ments to  duplicate  the  response  of  a flawed  structural  ele- 
ment to  fatigue  loading.  It  is  desired  to  have  a procedure 
which  does  not  require  extensive  additional  testing  when- 
ever a new  laminate  layup  is  used.  The  procedure  utilized 
herein  is  based  on  the  use  of  nondestructive  inspection  (NDI) 
methodology  to  detect  and  measure, critical  flaws.  This  is 
supplemented  by  an  analysis  to  predict  the  effect  of  this 
flaw  size  on  the  structural  response  under  given  loading 
conditions.  This  report  describes  the  analytical  and  experi- 
mental studies  to  develop  such  a procedure  and  to  provide 
data  for  correlation  purposes. 

Many  types  of  flaws  can  exist  in  a composite  structural 
element.  In  order  to  develop  the  required  methodology,  a 
through-the-thickness  notch  and  an  interlaminar  disbond  were 
chosen  as  two  typical  flaws  whose  effect  on  the  residual 
strength  and  the  residual  lifetime  of  the  structural  element 
was  to  be  studied.  For  each  of  these  flaws,  the  following 
tasks  were  carried  out: 

1.  Exploring  NDI  techniques  to  detect  and  measure  pos- 
sible strength  and  service  life  limiting  flaws  in 
graphite  fiber  reinforced  epoxy  composites. 
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2.  Developing  an  analysis  in  combination  with  the 
supplemental  experimental  data  to  quantify,  the 
severity  of  the  flaw. 

3.  Conducting  representative  fatigue  tests  and  monitor 
ing  the  growth  of  the  flaw  using  NDI  techniques. 

4.  Incorporating  this  into  an  NDE  methodology  to  pre- 
dict quantitatively  the  residual  lifetime  and  the 
residual  strength  of  the  structural  element. 

5.  Carrying  out  an  experiment/analysis  correlation 
study  to  verify  the  predictions  made  by  the  NDE 
methodology. 
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SCOPE  OF  THE  PROBLEM 


Defects  in  fiber-reinforced  composite  laminated  struc- 
tural components  originate  during  fabrication  and  during  ser- 
vice. Fabrication  or  "birth  defects"  may  occur  during  the 
preparation  of  the  prepreg  (improper  yarn  spacing,  broken 
filaments) , laminate  layup  and  cure  (matrix  voids,  resin-rich 
regions,  interlaminar  damage),  and  structural  component  assem- 
bly (improper  drilling,  machining,  potential  flaw  growth 
regions  such  as  bolted  joints,  etc.).  Service  damage  will 
occur  due  to  a fatigue  load  spectrum  (which  is  random  in 
nature) , impact  or  FOD  (Foreign  Object  Damage) , and  environ- 
ment (UV,  moisture,  corrosion,  lightning,  etc.).  The  vari- 
ous types  of  flaws  that  may  originate  and  propagate  during 
the  service  life  of  a composite  structural  component  are 
listed  in  table  1.  These  result  in  a large  number  of  poten- 
tial failure  modes. 

The  initiation  and  growth  of  flaws  is  dependent  upon 
the  imposed  loads  and  the  laminate  construction.  The  vari- 
ous types  of  loads  that  affect  flaw  initiation  and  growth 
may  be  categorized  as:  inplane  loads  (tensile,  compressive, 
shear,  and  biaxial) ; transverse  loads;  and  thermal  loads. 

For  all  of  these  loads,  it  is  necessary  to  consider  their 
stochastic  nature.  Also,  it  is  necessary  to  consider  envi- 
ronmental effects.  The  laminate  characteristics  that  influ- 
ence the  flaw  growth  are  the  material  properties  of  the 
individual  laminae  (dependent  on  the  fiber  orientation) , 
the  stacking  sequence  of  the  layers,  the  geometry  of  the 
laminate,  and  the  location,  size,  and  shape  of  the  flaw. 

A knowledge  of  the  various  types  of  flaws  that  may  ini- 
tiate and  grow  in  a composite  structural  element  and  the  fac- 
tors that  influence  these  flaws  identifies  a problem  of 
enormous  complexity.  The  objective  of  the  present  program 
is  to  choose  representative  flaws  that  are  considered  crit- 
ical to  composite  structures  and  to  develop  a methodology 
by  which  the  damage  tolerance  of  the  flawed  structure  can 
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fcs  quantified.  The  damage  tolerance  of  a flawed  composite 
structural  element  will  be  quantified  here  in  terms  of  the 
residual  strength  and  the  residual  modulus  of  the  element 
when  subjected  to  known  service  conditions.  The  objective 
is  met  by  combining  an  analysis  of  the  flawed  composite 
structure  with  an  experimental  program  that  supplements  the 
analysis.  The  undesirable  alternative  of  an  extensive 
experimental  program  utilizing  spectrum  fatigue  loading 
would  make  the  cost  prohibitive. 

In  order  to  arrive  at  a definition  of  a critical  flaw, 
it  is  necessary  to  formulate  analytical  and/or  semi-erapirical 
models  which  are  capable  of  predicting  the  direction  and 
rate  of  growth  and  the  critical  size  (at  failure)  of  the  dif- 
ferent types  of  flaws  such  as  surface  flaws,  slit  notches 
through  the  thickness,  and  delaminations.  Due  to  the  com- 
plexity of  the  problem  of  failure  in  composite  laminates, 
it  is  impossible  for  a single  analytical  model  to  encompass 
the  whole  spectrum  of  fracture  and  fatigue.  A logical  ini- 
tial step  is,  therefore,  to  begin  with  simpler  flaw  geome- 
tries such  as  a through-the-thickness  slit  notch  and  an 
interlaminar  disbond.  Based  on  the  analytical  and  experi- 
mental correlation  studies  and  the  NDE  tests  for  these  flaws, 
confidence  can  be  developed  for  similar  studies  for  other 
flaws. 

The  slit  notch  study  was  chosen  because  the  correspond- 
ing analysis  to  predict  the  failure  modes  and  loads  in  a 
notched  laminate  due  to  known  service  loading  conditions  was 
already  available  at  the  beginning  of  the  program.  The 
analysis  for  a composite  structure  with  a delamination  at 
arbitrary  locations  was  developed  under  the  present  program 
and  the  proposed  NDE  methodology  was  demonstrated  for  this 
defect  in  a qualitative  manner. 

The  problem  of  "critical  flaw  evaluation"  can  be  sum- 
marized as: 
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(i)  Determining  what  (type  of  flaw) , where  (location 
of  the  flaw),  when  (inspection  interval),  and  how 
(type  of  NDI  technique)  to  look  for  flaws. 

(ii)  Formulating  mathematical  models,  for  known  loading 
conditions,  to  predict  in  situ  critical  flaw  sizes,, 
modes  of  failure,  residual  lifetime,  and  residual 
strength  for  the  chosen  types  of  flaw.  This  is 
done  with  the  help  of  supplemental  experimental 
data. 

(iii)  Using  NDI  techniques  to  measure  the  actual  growth 
of  the  flaw  under  service  conditions. 

(iv)  Incorporating  (ii)  and  (iii)  into  an  NDE  methodology 
to  quantify  the  in  situ  residual  strength  and  resid- 
ual lifetime. 

(v)  Carrying  out  a correlatior  study  on  the  NDE  predic- 
tions. 

(vi)  Exploring  alternative  measurement  techniques  to 
improve  the  methodology. 
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ANALYTICAL  METHODOLOGY 

The  understanding  of  the  damage  tolerance  of  composite 
structures  is  a complex  problem.  A possible  method  of 
approach  to  this  problem  is  to  take  the  design  and  analysis 
methods  developed  for  metals  and  carry  them  over  directly 
for  composites.  If  the  method  works,  an  advance  in  the 
state  of  the  art  has  been  obtained  at  minimum  effort.  Thus, 
for  both  the  fatigue  and  fracture  problem  of  composites,  it 
is  reasonable  to  start  by  attempting  to  utilize  the  methods 
of  fracture  mechanics  developed  for  metals.  However,  one 
should  be  prepared  to  accept  the  fact  that  this  translation 
of  art  may  not  always  be  successful.  This  is  because  the 
fracture  and  fatigue  processes  in  fiber-reinforced  laminates 
are  expected  to  be  functions  of  some  of  their  unique  charac- 
teristics. These  include  heterogeneity,  statistical  varia- 
bility, microstructural  crack  growth,  and  the  three- 
dimensional  nature  of  the  stress  field. 

Several  approaches  to  the  predictions  of  fracture  or 
crack  growth  in  composite  materials  may  be  taken;,  namely,  a 
classical  fracture  mechanics  (CFM)  approach  applied  at  a 
micromechanical  level,  a classical  fracture  mechanics  approach 
applied  to  an  effective  material  on  a macromechanical  level, 
or  the  quasi-heterogeneous  minimechanics  approach  utilized 
in  the  earlier  studies  performed  at  Materials  Sciences 
Corporation  (e.g.  ref.  1). 

Classical  fracture  mechanics  applied  on  a micromechanical 
level  is  capable  of  treating  a composite  as  a heterogeneous 
medium  and  appears  in  concept  to  be  a valid  approach.  How- 
ever, the  extreme  complexity  of  the  required  analysis  places 
limitations  on  its  practicality.  Classical  fracture  mechanics 
applied  on  a macroscopic  scale  is  useful  for  certain  types  of 
problems.  The  methodology  of  CFM  may  be  valid  when  the  crack 
propagation  is  self-similar.  Thus,  the  scope  of  this  approach 
is  limited  because,  in  many  laminates,  planes  of  weakness, 
rather  than  the  initial  flaw  geometry,  determine  the  direction 
of  crack  propagation. 
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A slit  notch  and  an  interlaminar  delamination  in  a 
composite  laminate  have  been  identified  as  two  typical  flaws 
of  interest.  A static  failure  model  for  a laminate  with  a 
slit  notch  has  already  been  developed  at  Materials  Sciences 
Corporation  (refs.  1 to  5) . The  model  for  a delaminated 
composite  beam  was  developed  under  the  present  program.  These 
two  static  failure  models  are  discussed  below: 

STATIC  FAILURE  MODEL  FOR  A LAMINATE  WITH  A SLIT  NOTCH 

A brief  summary  of  the  salient  features  of  the  existent 
static  failure  model  for  a notched  composite  laminate  is  pre- 
sented in  this  section.  The  behavior  of  notched  laminates 
exhibits  some  features  which  have  not  been  observed  in  metals. 
They  are: 

(i)  The  direction  of  crack  growth  is  not  always  collinear 
with  the  initial  crack  and,  for  many  laminates,  it 
is  a function  of  notch  size,  environment,  and  stack- 
ing sequence. 

(ii)  The  measured  strength  reduction  resulting  from  the 
introduction  of  a circular  hole  in  a laminate  does 
not  correlate  with  that  obtained  from  the  stress 
concentration  factor  determined  from  anisotropic 
plate  theory.  Also,  the  reduction  in  strength  in- 
creases with  hole  radius;  that  is,  stress  concentra- 
tion factor  increases  with  increase  in  hole  size  in 
an  infinite  plate.  This,  too,  is  contrary  to  the 
results  obtained  for  a homogeneous  anisotropic  mate- 
rial. 

(iii)  In  contrast  to  conventional  metals  where  notch  sen- 
sitivity increases  as  the  unnotched  strength  in- 
creases, certain  fiber  composites  exhibit  increased 
notch  toughness  with  increasing  unnotched  strength. 

(iv)  Some  notched  laminates  show  a lack  of  sensitivity 
to  the  size  and  shape  of  notches. 
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(v)  The  residual  static  strength  of  some  notched 

laminates  subjected  to  fatigue  loading  is  equal 
to  or  greater  than  their  static  strength.  This 
is  again  in  contrast  to  the  high-strength  metal 
alloys  wherein  a fatigue  loading  increases  the 
growth  of  cracks,  leading  to  unstable  fracture. 

The  basic  static  failure  model  for  axial  and  transverse 
failure  of  notched  composites  was  developed  in  references  1,  2, 
3,  and  4.  The  model  for  the  static  failure  of  a notched  com- 
posite laminate  (ref.  1)  is  illustrated  in  figure  1.  The 
laminate  is  assumed  to  be  under  a tensile  stress  in  the  x 
direction.  A notch  of  width  "a"  is  centered  in  the  specimen 
at  x = 0.  The  central  core  region  tends  to  "pull  out"  from 
the  notch  area  due  to  the  applied  tensile  loading.  The  "pull 
out,"  however,  is  restrained  by  shear  stresses  between  the 
core  region  and  the  adjacent  intact  material.  These  shear 
stresses  generally  result  in  a region  of  high  shear  strain 
parallel  to  the  loading  direction.  Immediately  adjacent  to 
the  notch  core,  an  overstressed  region  of  average  stress  con- 
centration of  width  "a0u  is  assumed  to  exist.  Everywhere  out- 
side the  core  and  overs tressed  regions,  the  laminate  is  uni- 
formly strained.  Shear  strain  due  to  core  pullout  is  assumed 
to  extend  in  the  "y"  direction  over  a region  three  times  the 
size  of  the  overstressed  fiber  region  "ag" . This  assumption 
is  based  on  the  premise  that  steep  shear  stress  gradients  in 
the  core  and  the  average  material  adjacent  to  the  overstressed 
region  would  exist  over  a characteristic  dimension  equal  to 
"a0".  The  laminate  axial  shear  stress  - shear  strain  curve  is 
assumed  to  be  linear  elastic  - perfectly  plastic. 

The  static  failure  analysis  has  the  capability  of  com- 
puting the  axial  inelastic  length,  the  axial  crack  length, 
and  the  maximum  overstress  in  the  material  adjacent  to  the 
notch  for  a given  applied  laminate  tensile  stress.  It  also 
monitors  the  transverse  crack  propagation  mode.  The  analysis 
can  also  predict  off-axis  cracking  in  an  approximate  fashion 
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(see  ref.  4).  In  addition  to  axial,  transverse,  and  off-axis 
cracking,  interlaminar  delamination  may  also  originate  around 
the  notch  precipitating  laminate  failure.  Modifications  to 
the  failure  model  to  include  interlaminar  effects  were  incor- 
porated in  reference  5.  The  modified  analysis  can  predict 
axial  damage  at  the  lamina  level. 

Given  the  laminate  mechanical  properties  and  failure 
stresses  and  strains,  the  unknown  parameter  in  the  analysis 
is  "a0"  which  defines  the  extent  of  the  average  stress  concen- 
tration region.  Other  investigators  have  also  recognized 
the  need  to  define  the  dimension  of  this  region.  Examples  are 
the  "intense  energy  region"  (ref.  6)  and  the  distances  asso- 
ciated with  the  "point  stress"  and  "average  stress"  criteria 
(ref.  7).  Accurate  determination  of  "a0"  is  an  important 
aspect  of  the  predictive  capability  of  the  static  fatigue 
failure  model.  An  empirical  form  (based  on  existing  experi- 
mental data)  for  predicting  "ag"  as  a function  of  the  notch 
size  and  laminate  construction  has  been  postulated  (ref.  8). 
The  extent  of  the  transverse  damage  region  "ag"  is  postulated, 
based  on  some  physical  considerations,  while  the  lamina  axial 
damage  regions  are  determined  from  the  analysis  (ref.  5). 

The  existent  static  failure  analysis  for  notched  com- 
posite laminates,  described  above,  will  be  used  to  predict 
the  failure  modes  and  strengths  of  the  notched  laminates 
tested  experimentally  under  this  program. 

STATIC  FAILURE  MODEL  FOR  A LAMINATED  BEAM  WITH  AN  INTER- 
LAMINAR DELAMINATION 

Classical  fracture  mechanics  is  employed  to  develop  a 
static  failure  model  for  a tip-loaded  debonded  cantilever 
beam.  The  methodology  deployed  herein  has  been  used  by 
others  in  solving  similar  problems  (refs.  9 to  15).  A few 
assumptions  are  made  to  reduce  the  problem  to  a tractable 
form  and  to  gain  an  initial  knowledge  of  the  behavior  of 
delaminated  beams.  The  most  important  among  these  is  the 


-10- 


NADC-76228-30 


assumption  that  the  debond  propagates  along  the  interface 
in  a collinear  (mode  II)  fashion.  This  assumption  is  justi- 
fied by  testing  a tip-loaded  cantilever  beam  with  an  implanted 
delamination. 

The  static  failure  model  analyzes  the  debonded  beam  as 
four  separate  beam  elements  joined  together  at  the  crack  tips 
with  the  proper  boundary  and  matching  conditions  (fig.  2) . 

Each  beam  element  is  treated  as  a Timoshenko  beam  to  incor- 
porate shear  deformation  effects.  The  corresponding  equa- 
tions, boundary  conditions,  and  the  displacement  solutions 
for  the  four  elements  are  presented  in  Appendix  A.  The  total 
strain  energy  in  the  tip- loaded,  cantilevered,  debonded  beam 
is  obtained  using  Clapeyron's  theorem.  The  strain  energy  is 
dependent  on  the  location  and  the  length  of  the  delaraination, 
the  magnitude  of  the  applied  load,  and  the  material  and 
geometric  properties  of  the  laminated  beam  (Appendix  A) . 

The  criticality  of  the  delamination  is  determined  through 
Griffith's  energy  balance  criterion  (ref.  16).  This  classical 
fracture  mechanics  approach  states  that  the  loss  in  the  total 
strain  energy  due  to  an  incremental  change  in  the  crack 
length  is  equal  to  the  surface  fracture  energy  that  is  neces- 
sary to  create  the  corresponding  surface  area.  Inherent  in 
the  application  of  this  criterion  is  the  knowledge  of  the 
measure  of  ya#  the  specific  adhesive  surface  fracture  energy 
required  to  create  unit  surface  area  through  the  propagation 
of  the  delamination.  The  change  in  strain  energy  is  computed 
by  evaluating  the  energy  for  two  different  crack  lengths 
using  the  above  model. 

The  application  of  Griffith's  energy  balance  criterion 
leads  to  a closed-form  solution  for  P„  /fY-  as  a function  of 
the  beam  properties  (Appendix  A) . The  solutions  for  two  loca- 
tions of  the  debond  in  a beam  having  the  laminate  configura- 
tion used  in  the  tests  are  presented  in  figure  3.  B is  the 
magnitude  of  the  tip  load  at  which  propagation  of  the  delamina- 
tion of  size  a/L  is  imminent.  An  experimental  data  point, 
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used  with  this  solution,  will  suffice  to  estimate  the  magni- 
tude of  y . 

a 

It  is  shown  in  Appendix  A that  the  critical  value  of 
the  tip  load  is  a function  of  the  location  and  length  of 
delamination,  the  material  and  geometric  properties  of  the 
beam,  and  the  specific  adhesive  surface  fracture  energy  for 
the  given  adherend-adhesive  combination.  The  propagation  of 
the  delamination,  along  the  interface,  may  take  place  in  any 
of  the  following  directions  depending  on  the  material  and 
geometric  data: 

(i)  propagation  of  the  delamination  into  element  1; 

(ii)  propagation  of  the  delamination  into  element  2; 

(iii)  propagation  of  the  delamination  into  elements 

1 and  2 simultaneously. 

The  critical  load  for  each  case  is  obtained  and  the  direction 

in  which  the  delamination  propagates  is  governed  by  the  least 

of  the  three  P values, 
cr 


Finite  Element  Analysis  of  a Delaminated  Beam 

The  design  of  beam  specimens  to  study  the  growth  of 
delaminations  under  transverse  shear  conditions  introduces 
a few  problems  that  cannot  be  handled  by  beam  theory.  One 
is  the  shear  stress  concentration  directly  under  the  load  and 
the  resultant  change  in  the  shear  stress  distribution  from 
that  predicted  by  Euler  Bernoulli  beam  theory  (fig.  4).  This 
localized  stress  redistribution  would  lead  to  a lower  shear 
stress  at  the  crack  tip  if  the  debond  is  located  at  the  mid- 
plane or  in  the  lower  half  of  the  beam.  Consequently,  a 
larger  load  would  have  to  be  applied  to  exceed  the  shear 
strength  of  the  debonded  beam  and  make  the  delaminai  ion  propa- 
gate. This  could  cause  premature  failure  in  the  flexural  mode 
in  specimens  designed  to  measure  criticality  of  delaminations. 

Static  tests  on  t04/±452^2s  beam  specimens  were  conduc- 
ted under  this  program  to  measure  the  static  failure  loads. 
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with  and  \ ithout  interlaminar  delaminations,  and  to  obtain 
a quantitative  measure  of  the  specific  adhesive  surface  frac- 
ture energy.  The  results  are  listed  in  table  2.  It  can  be 
seen  that  the  debond  propagated  (in  an  unstable  manner)  only 
in  one  of  the  tests.  All  the  other  specimens  failed  in 
flexure,  perhaps  due  to  the  effect  shown  in  figure  4.  This 
prompted  the  use  of  a finite  element  analysis  to  design  a 
delaminated  beam  such  that  debond  propagation  precedes  flex- 
ural failure.  The  ANSYS  code  was  used  for  this  purpose. 

A plane  strain  finite  element  analysis  of  a debonded 
beam  was  carried  out  to  ensure  that  the  specimen's  shear 
strength  is  exceeded  at  a load  level  where  the  flexural 
strength  is  not.  The  delaminated  beam  was  modeled  as  shown 
in  figure  5 and  6.  Input  and  output  characteristics  for 
the  constant  strain  elements  are  shown  in  figures  7 and  8. 

The  material  properties  of  the  elements  in  the  different 
layers  are  shown  in  table  3.  The  thickness  of  the  beam  was 
scaled  up  by  a factor  of  four  to  avoid  numerical  problems 
due  to  small  element  aspect  ratios. 

Figure  9 shows  the  transverse  shear  stress  variation 
along  the  interface  (midplane)  containing  a debond.  Figure 
10  shows  the  shear  stress  variation,  through  the  thickness, 
at  the  crack  tips.  For  the  debonded  beam  model  in  figure  5, 
it  is  assumed  that  the  maximum  bending  stress  (200  ksi  to 
break  the  fibers  in  the  0°  layer  at  the  top)  is  twenty  times 
larger  than  the  magnitude  of  the  shearing  stress  that  causes 
delamination  (about  10  ksi) . Beam  theory  and  results  from 
figures  9 and  10  s.^ow  that  the  tip  load  that  causes  delamina- 
tion is  about  one-fifth  of  the  value  that  causes  bending 
failure.  Hence  delamination  precedes,  and  possibly  precludes, 
bending  failure.  It  is  seen  from  both  the  figures  that  the 
shear  stress  at  the  crack  tip  nearer  the  applied  load  and 
farther  from  the  support  (x  = + a)  is  a maximum.  The 

debond  should,  therefore,  tend  to  propagate  toward  the  free 
end  of  the  tip-loaded,  cantilevered  beam.  A second  set  of 
experiments  were  conducted  on  l (04/+452/+452/04) S1 s beam 
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specimens  with  a reduced  test  section  length  in  order  to 
guarantee  propagation  of  the  debond  at  load  levels  well 
below  laminate  flexural  strength. 

FATIGUE  ANALYSIS 

The  most  direct  way  to  determine  laminate  tensile  and 
shear  behavior  under  cyclic  loads  is  to  perform  fatigue 
tests  on  unnotched  specimens.  This  method,  though  straight- 
forward, would  require  a large  number  of  tests  for  each 
laminate  arid  the  results  would  be  applicable  only  for  the 
specific  laminate  that  was  tested.  Each  new  laminate  would 
require  additional  tests.  Therefore,  a method  of  generating 
laminate  fatigue  behavior  from  lamina  fatigue  properties  is 
used. 


General  Philosophy 

On  fatigue  loading  of  a notched  or  debonded  laminate, 
high  axial  and  shear  stresses  near  the  crack  tip  cause  the 
material  in  the  vicinity  of  the  crack  to  be  degraded  much 
more  rapidly  than  the  rest  of  the  laminate.  This  spatial 
variation  of  the  property  degradation  in  the  laminate  is  a 
result  of  the  different  stress  states  in  the  individual 
laminae,  and  the  large  stress  concentration  at  the  crack  tips 
in  the  neighboring  laminae.  An  analytical  prediction  of 
the  three-dimensional  stress  state  in  a flawed  laminate  is 
a formidable  task,  and  thus  a simplified  numerical  solution 
is  used  as  the  alternative. 

The  parameters  that  influence  the  three-dimensional 
stress  state  in  a flawed  laminate,  and  hence,  the  property 
degradations,  are:  the  number  of  cycles  of  loading,  N; 
the  ratio,  S,  of  the  maximum  fatigue  load  to  the  static 
failure  load  of  the  virgin  (unnotched  and  without  delamina- 
tions) laminated  beam;  the  stress  ratio,  R,  chat  is  the  ratio 
of  the  minimum  to  the  maximum  fatigue  load;  the  frequency  of 
cyclic  loading,  <u;  and  the  ratio,  n , of  the  static  failure 

load  of  a flawed  beam  of  known  flaw  size  to  the  static 
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failure  load  of  the  virgin,  unflawed  laminated  beam.  For  a 
debonded  beam,  n is  a function  of  the  initial  debond  size, 
a/L,  and  can  be  obtained  from  the  static  failure  analysis 
(fig.  3). 

The  general  fatigue  philosophy  is  based  on  using 
unflawed  lamina  fatigue  data  to  predict  the  local  property 
changes  in  flawed  laminates  (refs.  1,  4,  5,  and  8).  Fatigue 
tests  are  to  be  conducted  on  unflawed  coupons  to  span  a 
wide  range  of  stress  states  and  the  desired  property  degrada- 
tions measured.  These  results  constitute  the  data  bank  for 
the  fatigue  philosophy.  A knowledge  of  the  complex  in  situ 
stress  state  in  each  lamina  can  then  be  used  in  conjunction 
with  the  experimental  data  to  predict  the  spatial  property 
degradation  in  each  ply  of  the  flawed  laminate.  NDT  can 
be  employed  to  measure  the  local  property  degradation  for 
correlation  studies. 

The  static  failure  analysis  for  a flawed  laminate  under 
a given  set  of  loading  conditions  requires  the  stiffness 
and  strength  properties  as  input  data.  Fatigue  failure 
analysis  is  carried  out  by  considering  a finite  number  of 
load  cycle  intervals,  and  using  the  degraded  properties  for 
each  interval  in  the  static  failure  analysis. 

The  fatigue  philosophy  outlined  above,  though  sound  in 
principle,  is  difficult  to  implement.  The  formidable  stress 
analysis  and  the  need  for  an  extensive  experimental  program 
to  collect  the  required  data  deem  a more  approximate  theory 
necessary.  An  empirical  approach,  used  in  the  present  pro- 
gram, is  discussed  below. 

Empirical  Approach  Used  in  the  Present  Program 

A simplified  fatigue  philosophy  (refs.  1,  4,  5,  and  8) 
is  used  in  the  present  program  to  mathematically  model  the 
property  degradation  in  an  empirical  manner.  The  spatial 
stress  variation  in  each  lamina  is  approximated  by  an  average 
stress  state  throughout  the  lamina.  This  approximates  the 
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spatial  variation  in  the  lamina  S value  (the  ratio  of  the 
maximum  fatigue  stress  to  the  static  strength  of  the  unflawed 
lamina)  by  an  average  S value  for  each  lamina  in  the  flawed 
laminate.  The  spatial  variation  of  the  property  degradations 
in  each  ply  is  thus  approximated  by  average  degraded  ply 
properties  to  simplify  the  analysis.  The  fatigue  degradation 
in  each  ply  in  a flawed  beam  can  then  be  assumed  to  be  a func- 
tion of  N and  the  laminate  S/n  value  which  is  the  ratio  of  the 
maximum  applied  fatigue  load  to  the  static  failure  load  of  the 
flawed  beam  of  known  initial  flaw  size. 

The  fatigue  analysis  is  carried  out  by  computing  the 
laminate  property  degradations  from  the  lamina  fatigue  data 
using  laminated  plate  theory.  This  is  done  for  the  four  beam 
elements  in  the  present  analysis.  The  lamina  fatigue  data  are 
obtained  from  experiments  on  unnotched,  perfectly  bonded  uni- 
directional specimens.  These  tests  are  carried  out  at  differ- 
ent S levels.  When  the  lamina  fatigue  data  are  used  in  the 
fatigue  failure  analysis  of  a flawed  beam,  the  n value  corres- 
ponding to  the  measured  flaw  size  is  computed  from  the  static 
failure  analysis  and  the  lamina  fatigue  data  corresponding  to 
the  laminate  S/n  value  are  used  in  the  laminate  analysis  to 
obtain  average  degradations  in  the  four  beam  elements.  Hence, 
for  different  initial  flaw  sizes,  different  levels  of  degrada- 
tions result  depending  on  the  ratio  of  the  maximum  fatigue  load 
to  the  static  failure  load  of  the  flawed  beam.  As  a further 
simplification  in  this  example,  the  rates  of  property  degrada- 
tion in  this  beam  are  assumed  to  be  constant  in  all  plies. 

Thus  the  effect  of  the  differences  in  property  degradation 
resulting  from  different  stress  levels  in  the  various  plies 
is  disregarded  to  simplify  the  example.  Different  rates  of 
degradation  are  used,  however,  for  the  different  moduli  of  the 
plies  (i.e.  E^,  E 2,  G12)*  These  rates  are  the  rates  of 
degradation  corresponding  to  the  initial  laminate  S/n  value. 

The  effect  of  the  change  in  laminate  S/n  on  the  degradation 
rates  is  ignored. 
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Figure  11  shows  the  variation  in  the  critical  load  with 
the  number  of  cycles  and  the  debond  size  for  a [ 04/± 4 52 3 4S 
beam  under  cyclic  loading.  The  debond  is  located  at  l^/ L 
**  1/3  and  it  is  at  the  laminate  middle  surface.  This  analysis 
does  not  have  the  capability  to  treat  stable  growth  of  the 
debond.  The  effects  of  crack  growth  will  be  discussed  sub- 
sequently. The  assumed  property  degradation  rates  are  stated 
on  the  figure.  These  rates  correspond  to  a laminate  S/n 
value  of  two-thirds  for  an  initial  debond  size  of  a/L  = 0.2. 

If  the  initial  debond  size  is  different,  different  curves 
(for  the  new  laminate  S/n  value)  result  for  N = 106  cycles 

7 

and  N = 10  cycles  in  figure  11. 

VIBRATION  ANALYSIS  OF  A DELAMINATED  BEAM 

The  degradation  in  the  material  properties  of  a flawed 
specimen  under  cyclic  loading  can  be  estimated  by  acoustic 
monitoring  or  ultrasonic  measurements.  Ultrasonic  C-scans 
were  used  extensively  in  the  present  study  to  detect  and 
measure  flaws.  Vibrational  monitoring  was  explored  as  a possi- 
ble NDT  for  future  use,  and  for  correlation  studies.  Fre- 
quency is  used  as  a parameter  to  study  the  criticality  of 
a flaw  (ref.  17),  because  the  reduction  in  the  stiffness 
as  well  as  the  growth  of  the  flaw  affect  the  natural  frequen- 
cies of  vibration  of  the  specimen.  An  analytical  methodology 
to  evaluate  the  reductions  in  the  frequencies  due  to  each 
cause  is  presented  in  Appendix  B for  a delaminated  beam. 

The  model  is  similar  to  the  static/fatigue  model  that  deter- 
mines the  criticality  of  interlaminar  delaminations  (Appendix 
A).  The  free  vibration  analysis  (see  fig.  12)  assumes  the 
beam  elements  to  be  Timoshenko  beams  that  include  shear 
deformation  and  rotary  inertia  effects. 

The  free  vibration  frequencies  for  the  debonded  beam  are 
obtained  by  computing  the  determinant  of  a 16x16  matrix  for 
many  assumed  values  of  frequencies,  and  numerically  inter- 
polating to  find  the  frequencies  corresponding  to  a null 
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determinant  (Appendix  B) . Normalizing  these  frequencies 
with  respect  to  the  corresponding  frequencies  of  vibration 
of  the  beam  without  any  delamination,  the  decrease  in  the 
frequencies  with  the  growth  of  debond  can  be  estimated.  Like- 
wise, for  any  size  of  the  debond,  the  reduction  in  the 
frequencies  with  cyclic  degradation  in  the  element  stiffnesses 
can  also  be  evaluated. 

Debonds  of  known  sizes  are  implanted  in  specimens  in  the 
laboratory  and  experimental  measurements  of  the  frequenices 
made  after  known  numbers  of  cycles.  Periodic  C-scans  of  these 
specimens  reveal  the  damage  growth  (stable)  with  cyclic  load- 
ing, making  feasible  a frequency -damage  size  relationship. 

The  experimentally  measured  frequencies  for  different  debond 
sizes  could  then  be  correlated  with  the  analytical  predic- 
tions. If  the  debond  grows  in  a stable  manner  till  the  onset 
of  catastrophic  propagation,  the  critical  state  could  be 
acoustically  monitored.  Acoustic  measurements  could  then  be 
used  as  a viable  tool  to  predict  the  criticality  of  a flaw 
(a  slit  notch,  delamination,  etc.)  in  laminated  structures. 

All  the  vibration  tests  were  conducted  on  [ (04/+452) S1 s 
beam  specimens.  The  first  four  natural  vibration  frequencies 
for  these  beams  were  obtained  with  no  interlaminar  defect 
(table  4) . The  same  frequencies  were  also  measured  for  beam 
specimens  with  1.27  cm  and  2.54  cm  delaminations  implanted 
in  them  (tables  5 and  6) . These  results  are  compared  with 
those  predicted  by  the  vibration  analysis  in  table  7.  While 
the  results  for  an  unflawed  beam  demonstrate  good  correlation, 
those  for  a debonded  beam  with  an  a/L  of  0.2  are  not  encour- 
aging and  show  reasonable  agreement  only  for  the  fourth  fun- 
damental frequency.  This  fact  calls  for  further  investiga- 
tion to  assess  whether  the  vibration  tests,  when  understood 
well,  will  serve  as  a useful  tool  for  NDE  of  damage  growth 
in  flawed  beams. 
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NDE  METHODOLOGY  TO  PREDICT  FLAW  CRITICALITY 

In  evaluating  the  integrity  of  structural  elements,  non- 
destructive inspection  focuses  upon  detection  of  flaws  and 
the  determination  of  their  sizes,  while  non-destructive  evalua- 
tion techniques  translate  these  into  a measure  of  the  cri- 
ticality of  the  flaws.  Allowable  limits  for  flaw  size  and 
damage  growth  rates  must  be  determined  to  establish  inspec- 
tion intervals  for  in-service  flawed  strvxctural  elements,  and 
to  determine  flaw  characteristics  that  make  repair  mandatory. 

To  this  end,  an  NDE  methodology  is  developed  herein.  It  is 
built  around:  (i)  a static/fatigue  failure  analysis;  (ii) 
experimental  results  on  flawed  specimen  response;  and  (iii) 
an  NDE  methodology  to  assess  the  criticality  of  a flaw  using 
NDT. 

The  analytical  results  predicted  by  the  static/fatigue 
failure  model  (figs.  3,  11)  form  the  basis  of  the  NDE  method- 
ology. For  a defect  of  a known  initial  size  and  location, 
a family  of  curves  describes  the  variation  of  the  critical 
load  with  the  defect  size  for  different  values  of  N (fig.  11) . 
As  the  number  of  cycles  (N)  of  known  loading  conditions 
increases,  the  spatial  degradation  of  the  material  properties 
causes  a lowering  of  the  critical  load  at  which  the  propa- 
gation of  a delamination  of  known  size  is  imminent  (fig.  11). 
The  crack  growth  is  stable  if  additional  load  is  required  to 
increase  the  crack  size,  and  unstable  if  an  increase  in  a/L 
results  in  a reduction  of  the  critical  load.  The  present 
model  for  an  interlaminar  defect  in  a beam  yields  values  of 
critical  load  at  which  the  imminent  propagation  of  the  debond 
is  unstable  by  nature.  Unstable  crack  growth  is  represented 
in  figures  3 and  11,  where  the  critical  load  decreases  with  a 
growth  in  the  delamination.  If  an  applied  load  reaches  the 
critical  value  for  the  existent  length  of  debond,  the  debond 
tends  to  grow  indefinitely.  Hence,  the  above-mentioned 
family  of  curves  describes  the  catastrophic  debord  growth 
for  each  value  of  N.  This  is  of  major  importance  to  the  NDE 
analyst. 
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For  a known  size  and  location  of  a debond  in  a delam- 
inated beam,  curves  of  "P  /J~y"  versus  a/L"  can  be  drawn  for 

cr  a 

various  values  of  N as  shown  in  figure  11.  As  discussed 
earlier,  these  curves  are  all  calculated  based  upon  the 
assumption  that  property  degradation  occurs  throughout  the 
laminate  at  constant  rates  (as  shown  in  figure  11) . In  actu- 
ality, the  degradation  rates  depend  upon  the  crack  size  and 
the  stress  ratio  S/n.  For  these  calculations,  the  degrada- 
tion rates  used  were  those  for  a/L  = 0.2  and  S/n  » 2/3. 

From  curves  of  this  type,  the  static  failure  load,  Pstat^c» 
for  any  initial  debond  size  can  be  obtained  from  the  N = 1 
curve.  For  illustrative  purposes,  pstatic/4  Ya  for  an  a^L 
value  of  0.2  is  selected  in  figure  11.  If  the  laminate  S 
level  for  fatigue  loading  is  chosen,  the  maximum  cyclic  load, 
^fatigue ' can  be  determined  as  Pfatique  for  the 

chosen  laminate  s value  corresponds  to  a particular  laminate 
S/n  value  for  the  known  debond  size.  pfatigUe  for  s/n  = 2^3 
is  shown  in  figure  11.  The  line  pcr/jY^  = pfatigue/jrYa 
intersects  the  family  of  failure  curves  at  various  points. 

The  points  of  intersection  of 'Pfati  e with  the  familY  of 
curves  can  be  plotted  as  an  "a/L  versus  log  N"  curve  (fig. 

13).  The  lower  solid  curve  in  figure  13  defines  the  life- 
time to  catastrophic  propagation  of  an  initial  debond  of  tjie 
specified  size.  The  upper  solid  curve  differs  from  this 
only  in  the  assumption  that  the  surface  energy  remains  con- 
stant while  the  stiffnesses  degrade. 

All  of  these  calculations  result  in  unstable  crack  propa- 
gation. Experimentally,  however,  a stable  crack  propagation 
is  observed.  Treatment  of  this  phenomenon  would  greatly 
complicate  the  analysis,  since,  strictly  speaking,  once  the 
debond  grows,  the  n value  (and  hence  the  laminate  S/n) 
changes.  Thus,  different  fatigue  data  are  required  and  a 
new  family  of  "p__/TyI  versus  a/L  for  various  N"  curves  should 
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be  used.  This  is  not  done  in  figure  13  for  reasons  explained 
in  the  section  "Empirical  Approach  used  in  the  Present  Pro- 
gram. " The  initial  debond  size  of  a/L  = 0.2  is  the  lows*  . 


NADC-76228-30 


bound  on  the  a/L  value.  (see  fig.  11)  is  the  number  of 

cycles  after  which  an  initial  debond  of  size  a/L  = 0.2  will 
propagate  catastrophically  without  any  intermediate  growth. 

Nf  is  a hypothetical  upper  bound  on  N,  since  in  reality,  the 
damage  may  grow  with  N,  thus  shortening  lifetime.  Point  B 
corresponds  to  that  debond  size  at  which  a static  load  of 
magnitude  equal  to  Pfatigue  precipitates  a catastrophic 
failure  (N  = 1) . On  fatigue  loading,  the  debond  in  a lam- 
inated beam  appears  to  grow  in  a stable  manner,  as  shown 
by  curve  AED  in  figure  13,  until  the  {a/L,  N)  combination 
corresponds  to  a catastrophic  situation.  The  use  of  the 
curves  of  figure  11  is  rigorously  correct  only  if  the  prop- 
erty degradation  rate  is  insensitive  to  debond  size.  Point  D 
in  figure  13  then  represents  a typical  (a/L,  N)  combination 
at  which  catastrophic  debond  growth  precipitates  after  stable 
crack  growth. 

Knowing  the  location  and  initial  size  of  the  debond, 
and  the  chosen  cyclic  maximum  load,  the  boundary  for  cata- 
strophic propagation  of  the  delaraination  is  obtained.  This 
boundary  is  drawn  on  an  "a/L  versus  N"  plot  for  every  value 
of  initial  debond  size.  The  procedure  for  this  follows  the 
steps  outlined  below: 

1.  Determine  the  initial  debond  size  (a/L  and  its 
location. 

2.  From  the  static  failure  analysis,  obtain  a "p  //v 

CX*  cl 

versus  a/L"  curve  for  the  debonded  beam. 

3.  Given  the  maximum  applied  cyclic  load  and  the 
static  failure  load  of  the  unflawed,  .virgin  beam, 
determine  the  laminate  S/n  value  for  the  initial 
debond  size  from  the  static  failure  analysis  results. 

4.  Using  lamina  degradation  data  corresponding  to  the 
initial  laminate  S/n  value,  obtain  a family  of 
curves  defining  catastrophic  growth  of  defect  for 
different  values  of  N. 
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5.  Corresponding  to  the  maximum  cyclic  load,  obtain  the 
critical  values  of  a/L  (greater  than  or  equal  to  the 
initial  value)  for  different  values  of  N. 

6.  Plot  these  points  on  an  "a/L  versus  N"  plot,  and 
draw  a curve  through  them  to  define  the  ^catastrophic 
debond  growth  boundary. 

Figure  13  presents  the  instability  boundary  for  a 

^°4/i4^2>/,+4^2^®4^  s laminated  beam  with  a debond  in  the  mid- 
plane. BFLGC  is  the  curve  that  defines  unstable  debond  growth 
or  failure.  While  the  outer  curve  assumes  no  degradation  in 
the  specific  adhesive  surface  energy,  the  inner  one  assumes 
a degradation  in  y,  that  is  of  the  same  order  as  the  shear 

a 

modulus  degradation,  based  on  the  reasoning  that  the  surface 
energy  is  affected  directly  by  a deterioration  in  the  shear 
resistance  at  the  adherend-adhesive  interface.  The  initial 
defect  size  (a/L  = 0.2)  is  labeled  A,  and  the  critical  defect 
size  for  static  failure  at  the  maximum  cyclic  load  is  labeled 
B.  Point  C denotes  the  hypothetical  number  of  cycles  the 
beam  can  go  through  until  failure  without  any  growth  of  the 
defect.  This  is  not  a physically  realizable  state  since  the 
damage  generally  grows  when  the  flawed  specimen  is  subject 
to  cyclic  loading. 

The  actual  growth  of  the  damage,  which  can  be  monitored 
ultrasonically , may  follow  a path  similar  to  that  shown  by 
the  discontinuous  line  AED.  From  an  NOE  point  of  view,  if 
the  specimen  is  inspected  after  a few  cycles  and  the  defect 
has  grown  to  the  extent  shown  by  E,  the  residual  strength  and 
lifetime  of  the  specimen  have  to  be  determined  to  decide  on 
the  need  for  mandatory  repair.  If  the  initial  debond  (a/L  -■ 
0.2)  has  grown  to  state  E,  after  about  5.75  x 106  cycles,  FOG 
becomes  the  new  envelope  that  bounds  stable  crack  growth. 

At  E,  if  the  fatigued  beam  specimen  were  tested  statically, 
from  the  measure  of  EF  and  the  static  failure  analysis  (fig. 

3) , the  residual  strength  in  the  beam  may  be  computed.  Also, 
an  overestimation  of  the  residual  lifetime  of  the  beam  may 


E to  remain  constant  when  cyclic  loading  on  the  delaminated 
beam  is  assumed  to  be  applied  for  an  additional  number  of 


cycles  (given  by  EG) . Hence,  from  figure  13,  it  is  seen  that 
EF  is  a measure  of  the  residual  strength  in  the  specimen  and 
EG  is  an  overestimation  of  its  residual  lifetime.  These  could 


then  be  compared  with  a quantified  criterion  for  mandatory 
repair  and  serve  as  an  assessment  of  the  criticality  of  the 
flaw. 
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EXPERIMENTAL  PROGRAM 

The  primary  objectives  of  the  experimental  program  are: 

(i)  to  develop  material  property  data  required  to  support 
the  analysis;  (ii)  to  produce  realistic  birth  and  service 
defects  in  the  composite  material  laminates  for  subsequent 
periodic  nondestructive  evaluation;  and  (iii)  to  conduct 
experiments  required  to  evaluate  the  accuracy  and  limitations 
of  the  developed  analytical  models  in  assessing  flaw  criti- 
cality in  composite  material  systems.  The  two  generic  de- 
fects to  be  examined  are  the  slit  notch  (through  the  thick- 
ness) and  the  interlaminar  debond  or  delamination. 

The  above-mentioned  objectives  are  met  through  the  fol- 
lowing experimental  tasks : 

(i)  Fabrication  of  AS/3501  Gr/epoxy  (flawed  and  un- 
flawed) laminates  and  conduct  ultrasonic  ND 
testing  of  laminates  after  fabrication. 

(ii)  Determination  of  static  (longitudinal  and  trans- 
verse tension,  inplane  and  interlaminar  shear) 
and  fatigue  (S-N  and  residual  strength/stiffness-N 
for  longitudinal  tension  and  inplane  shear)  proper- 
ties for  the  AS/3501  material  system.  These  data 
are  utilized  in  the  analysis  presented  earlier. 

(iii)  Performance  of  unnotched  and  notched  (slit  notch) 
laminate  static  and  fatigue  tests  to  determine 
damage  growth,  residual  strength,  and  lifetime. 

(iv)  Performance  of  unflawed  and  flawed  laminated 
beam  static/fatigue  tests  to  determine  load- 
deflection  characteristics  and  growth  of  debond. 

The  ND  tests  for  damage  detection  and  damage 
growth  determination  are  performed. 

(v)  ND  vibration  tests  on  the  laminated  beam  to  assess 
the  effect  of  debond  growth  on  the  first  four  fun- 
damental frequencies  of  vibrations. 
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SPECIMEN  FABRICATION 

All  test  specimens  were  fabricated  from  a single  graphite- 
epoxy  prepreg  system  (Hercules  AS/3501)  purchased  in  accord- 
ance with  the  acceptance  standards  of  the  McDonnell-Douglas 
Corporation.  All  graphite-epoxy  laminates  were  fabricated 
in  a controlled  environment  using  hand  layup  procedures  and 
processed  in  the  University  of  Delaware  autoclave  system, 
using  conditions  wh^ch  yield  properties  typical  of  the  given 
composite  system  (AS/3501) . Sufficient  acceptance  testing 
was  performed  to  assure  the  quality  of  composite  laminates 
fabricated.  After  fabrication,  samples  of  all  laminates 
were  supplied  to  the  NADC  technical  monitor  for  inspection 
(see  figs.  14 , 15). 

MECHANICAL  CHARACTERIZATION 

This  section  describes  the  mechanical  characterization 
test  program  required  to  characterize  the  Hercules  AS/3501 
composite  material  system. 


Fiber  Direction  Tension  Properties 


The  axial  tension  tests  were  conducted  on  specimens  of 
geometry  shown  in  figure  15.  The  test  specimen  geometry  was 
22.86  cm  (9.0  inches)  in  length  and  1.27  cm  (0.5  inch)  in 
width  with  a test  section  of  15.24  cm  (6.0  inches).  The 
failed  eight-ply  specimens  are  shown  in  figure  16  where  con- 
siderable axial  splitting  is  apparent.  The  average  mech- 
anical properties  determined  from  the  tests  (fig.  17  and 
table  8)  are  as  follows:  ultimate  tensile  strength  of  1703 
MPa  (247  ksi) ; Young's  modulus  (E^)  of  125  GPa  (18.2  x 106 
psi);  and  Poisson's  ratio  (v12)  of  0.28. 
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Transverse  Tension  Properties 

Test  specimen  geometry  for  the  transverse  tension  test 
is  also  shown  in  figure  15.  The  specimen  is  22.86  cm  (9.0 
inches)  in  length  and  2.54  cm  (1.0  inch)  in  width.  The 
failed  eight-ply  specimens  are  shown  in  figure  18.  Mech- 
anical properties  determined  from  the  transverse  tension 
tests  (fig.  19  and  table  8)  are:  a Young's  modulus  (E0) 
of  9.93  GPa  (1.44  x 10  psi)  and  an  ultimate  transverse  ten- 
sile strength  of  58.8  MPa  (8.53  ksi)  . 

Inplane  Shear  Properties 

The  inplane  shear  properties  for  the  graphite-epoxy  com- 
posite were  determined  with  the  (±45)  , tensile  coupon  test 
specimens  shown  in  figure  15.  Test  specimen  dimensions  were 
identical  to  those  of  the  transverse  tension  test  specimens 
described  in  the  previous  paragraph.  Inplane  shear  proper- 
ties determined  (fig.  20  and  table  8)  include  an  ultimate 
shear  strength  of  88.2  MPa  (12.7  ksi)  and  a shear  modulus  of 
5.45  GPa  (0.79  x 106  psi).  The  failed  eight-ply  [±45] 

s 

laminate  test  specimens  are  shown  in  figure  21. 

Interlaminar  Shear  Properties 

The  interlaminar  shear  properties  were  determined  (fig. 
22  and  table  9)  through  the  thick  laminate,  short  beam, shear 
tests.  The  test  specimens  were  beams  of  aspect  ratio  4/1 
carrying  a concentrated  load  at  midspan.  Shear  strain  was 
monitored  by  means  of  a strain  rosette  at  the  beam  quarter- 
span.  Specimen  geometry  was  1.27  cm  (0.5  inch)  in  thick- 
ness, 7.62  cm  (3.0  inches)  in  length,  and  0.76  cm  (0.3  inch) 
in  width.  The  short  beam  shear  specimens  are  shown  in  figure 
23.  The  test  apparatus,  with  the  beam  mounted  in  place,  is 
shown  in  figure  24. 
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STATIC  LAMINATE  TESTS 

Laminates  with  and  without  defects  have  been  studied 
experimentally. 


Unnotched  Laminate  Properties 

The  response  of  the  unnotched  [45/0/-45/0]  laminate 

s 

was  determined  for  tensile  loadings  of  the  laminate  specimen 
shown  in  figure  15.  The  specimen  geometry  was  30.4  cm 
(12.0  inches)  in  length  and  3.81  cm  (1.5  inches)  in  width, 
with  a t st  section  of  22.9  cm  (9.0  inches).  The  measured 
laminate  properties  (fig.  25  and  table  8)  included  an  ulti- 
mate tensile  strength  of  866.0  MPa  (125.7  ksi),  Young's 
modulus  of  72.6  GPa  (10.5  x 10^  psi) , and  Poisson's  ratio 
of  0.70.  The  failed  laminate  specimens  are  shown  in 
figure  26. 


Notched  Laminate  Properties 

The  response  of  the  (45/0/-45/0]  laminate  containing 

s 

a slit  notch  defect  was  investigated  next.  The  slit  notch 
was  ultrasonically  machined  into  the  specimen  (fig.  27).  The 

slit  notch  was  0.952  cm  (0.375  inch)  in  length  with  a notch 

t —2  —"5 

| txp  radius  of  1.9  x 10  cm  (7.5  x 10  inches).  The  static 

test  results  are  shown  in  figure  28.  The  average  ultimate 
strength  was  determined  (table  8)  to  be  292.4  MPa  (42.4  ksi) 
and  the  measured  Young's  modulus  was  72.4  GPa  (10.5  x 106 
psi) . The  failed  notched  laminates  are  shown  in  figure  29. 

Laminate  Flexural  Behavior 

In  order  to  assess  the  influence  of  an  interlaminar 
defect  upon  the  flexural  strength  of  a composite  laminate, 
defects  were  induced  by  a tubular  film  embedded  between  plies 
32  and  33  in  a 64-ply  (04/±452l4s  ^•am^nate*  Defect  lengths 
of  2.54  cm  and  1.27  cm  (1.0  inch  and  0.5  inch)  were  con- 
sidered. The  beam  was  2.54  cm  (1.0  inch)  wide  with  a span 


1 


t 

> 


A. 


NADC-7 6228-30 

of  7.62  cm  (3.0  inches).  A typical  interlaminar  defect  is 
shown  at  400X  magnification  in  figure  30.  Failure  resulted 
through  the  unstable  propagation  of  the  interlaminar  defect 
to  complete  failure.  Typical  failure  specimens  are  shown  in 
figure  31.  Average  failure  load  (table  10)  for  the  2.54  cm 
defect  was  2870  (292  kg)  and  4652  (466  kg)  for  the  1.27  cm 
defect.  Load-deflection  curves  for  an  unflawed  specimen  and 
one  with  an  implanted  debond  are  presented  in  figures  32  and 
33.  These  can  be  used  to  measure  the  bending  stiffness  of 
the  beam  specimens  for  correlation  studies. 

Birth  Defect  Detection 

The  ultrasonic  "C"  scan  technique  was  employed  to  examine 
each  laminate  after  fabrication  to  locate  birth  defects  intro- 
duced during  the  fabrication  process.  Figure  34  shows  a typ- 
ical "C"  scan  of  a [45/Q/-45/0]  laminated  panel.  Note  that 
fiber  orientations  and  regions  of  thickness  variation  are 
readily  observed.  Specimens  were  machined  out  of  this  panel 
to  be  used  in  the  various  tests.  To  account  for  damages  due 
to  the  machining  process,  the  ultrasonic  technique  was  again 
deployed,  after  each  specimen  was  fabricated  and  prior  to 
installation  of  instrumentation.  This  was  to  assess  total 
specimen  fabrication  damage  and  to  provide  a reference  stand- 
ard. Figure  35  shows  results  for  several  unnotched  laminate 
specimens. 

FATIGUE  AND  VIBRATION  TESTS 

Fatigue  tests  were  performed  both  on  unidirectional  ma- 
terial coupons  and  on  laminates.  Vibration  tests  were  con- 
ducted on  beams  containing  interlaminar  disbonds. 

Fiber  Direction  Tension  Fatigue  Behavior 

Axial  tension  specimens  were  Subjected  to  fatigue  load- 
ings of  R * 0.1  at  30  Hertz,  and  at  stress  levels  of  S * 0.67 
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and  0.80  (ratios  of  cyclic  maximum  to  static  failure  stress). 
All  specimens  subjected  to  80  percent  loadings  exhibited 
fatigue  failures  (table  11) . The  specimens  subjected  to  67 
percent  loadings  and  cycled  to  500,000  and  1,000,000  cycles 
exhibited  no  reduction  in  strength.  In  addition,  no  damage 
was  detected  by  ultrasonic  inspection. 

Inplane  Shear  Fatigue  Behavior 

The  inplane  shear  fatigue  damage  was  examined  by  subjec- 
ting [±45]  coupons  to  tension-tension  fatigue  loadings  for 
S = 0.45  and  0.55  at  30  Hertz,  and  R = 0.1.  Figure  36  shows 
the  difference  between  ultrasonic  traces  for  samples  before 
and  after  a fatigue  loading  of  S = 0.45  and  N = 10,000  cycles. 
The  figure  shows  the  free  edge  damage  in  the  specimen  due  to 
fatigue  loading.  The  corresponding  results  are  presented  in  ' 
table  12. 


Laminate  Tension-Tension  Fatigue  Behavior 

The  unnotched  [45/0/ -45/01  laminate  was  subjected  to 

s 

fatigue  loadings  of  S = 0.67  and  0.80  at  30  Hertz,  and  R = 0.1. 
The  results  are  presented  in  tables  13  and  14.  The  accumula- 
tion of  edge  damage  after  500,000  and  1,000,000  cycles  at 
S * 0.67  is  shown  in  figure  37.  The  figure  shows  that  signif- 
icant edge  damage  results  from  tensile  fatigue  loading  of  the 
specimen  and  that  ultrasonic  techniques  are  sufficient  for  its 
detection.  The  variation  in  the  experimentally  measured 
residual  strength  and  the  residual  modulus  of  the  unnotched 
[45/0/-45/0]s  laminate  with  the  number  of  cycles  of  fatigue 
loading  is  shown  in  figure  38. 

Ultrasonic  inspection  results  for  pristine  notched  lam- 
inate specimens  are  shown  in  figure  39.  After  these  specimens 

C g 

were  subjected  to  10  , 5x10  , and  10  cycles  (table  15), 
ultrasonic  inspection  was  carried  out  to  examine  the  damage 
in  the  vicinity  of  the  notch.  These  results  are  shown  in  fig- 
ures 40  and  41,  where  the  ultras^bic  frequencies  are  4.2  MHz  and 
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7.2  MHz , respectively.  Subsurface  damage  at  the  tip  of  the 
slit  notch  for  each  specimen  is  particularly  noticeable  in 
figure  41.  In  addition,  the  damage  area  is  seen  to  increase 
with  an  increase  in  the  number  of  cycles.  A magnified  view 
of  a fatigued  specimen  (X212B)  in  figure  40  is  shown  in 
figure  42  for  a better  picture  of  the  damage.  The  results 
in  table  15  are  plotted  in  figures  43  and  44  to  show  the  varia- 
tion in  the  residual  strength  and  the  residual  modulus  at  two 
S values  for  a notched  I45/0/-45/0]  laminate. 

o 

Laminate  Flexural  Fatigue  Behavior 

Graphite/epoxy  laminates  containing  implanted  interlaminar 
defects  were  subjected  to  flexural  fatigue  loadings  to  examine 
cyclic  propagation  of  interlaminar  defects.  The  flexural  ap- 
paratus is  shown  in  figure  45 , where  it  is  seen  that  the  load 
is  introduced  at  midspan  and  resisted  by  simple  supports. 

The  implanted  defect  is  at  quarter-span.  The  complete  fatigue 
set-up  is  shown  in  figure  46.  Typical  results,  presented  in 
figure  47,  show  that  the  interlaminar  defect  propagates  from 
the  point  of  intersection  of  the  loading  edge  of  the  defect 
and  the  free  edge  of  the  specimen.  In  addition,  these  results 
show  that  the  defect  propagates  parallel  to  the  45°  fibers  in 
the  laminate  although  the  defect  is  positioned  between  two  0° 
layers  in  the  laminate.  At  an  S level  of  0.5,  the  damage  is 
shown  to  grow  exponentially,  and  catastrophic  failure  occurs 
at  N - 90,000  cycles.  Again,  the  ultrasonic  technique  is 
shown  to  be  adequate  for  detection  of  the  damage  and  its  sub- 
sequent propagation.  The  damage  growth  measured  ultrasonically 
is  plotted  in  figures  48  and  49. 

Laminate  Flexural  Vibration  Characteristics 

The  (0,4/±452l 2a  laminated  beam  specimens  were  subjected 
to  vibration  tests  after  known  cycles  of  fatigue  loadings  to 
assess  the  influence  of  interlaminar,  damage  upon  the  material 
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frequencies.  The  vibration  apparatus,  with  the  specimen 
exhibiting  the  third  and  fourth  mode  shapes,  is  shown  in 
figures  50  and  51,  respectively.  Relatively  insignificant 
dependence  of  natural  frequency  on  interlaminar,  damage  was 
observed  in  the  experiments  (tables  4,  5 and  6). 
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DISCUSSION  OF  RESULTS 


SLIT  NOTCH  STUDY 

A slit  notch,  through  the  thickness  of  the  laminate, 
was  chosen  to  be  one  of  the  two  flaws  to  be  studied  under 
the  present  program  because  of  the  importance  of  through 
holes  and  of  the  availability  of  an  existing  analysis  to 
study  the  effect  of  a slit  notch  under  tensile  fatigue  load- 
ing. Notched  [45/0/-45/0]  AS/3501  laminates  were  chosen 

s 

and  the  virgin  specimens  C- scanned  (fig.  39) . Some  of  these 
specimens  were  used  to  characterize  the  notched  laminate 
properties  (figs.  28,  29  and  table  8).  The  others  were  used 
to  study  the  fatigue  behavior  of  the  notched  specimens 
(figs.  40.  41,  42  and  table  15). 

The  ratio  of  the  average  residual  Young's  modulus  after 

a known  number  of  cycles  to  the  static  (notched)  modulus  and 

the  ratio  of  the  average  residual  strength  to  the  static 

notched  strength  are  plotted  (figs.  43,  44)  for  two  S values. 

Allowing  for  a reasonable  scatter  in  the  experimental  results, 

the  Young's  modulus  seems  to  remain  almost  constant.  The 

residual  strength  increases  with  N for  an  S value  of  0.67 

(fig.  43).  But  for  S - 0.8,  the  residual  strength  increases 
5 

until  N = 5 x 10  cycles,  and  then  decreases. 

The  existing  static  failure  analysis  (ref.  5)  was  used 
to  predict  the  damage  growth  in  the  notched  [45/0/-45/01 

S 

AS/3501  specimen  under  a static  tensile  load.  The  analysis 
yields  a curve  relating  the  length  of  the  axial  damage  (in- 
cluding both  the  inelastic  zone  and  the  axial  crack)  to  the 
applied  load.  The  results,  shown  in  figure  52,  assume  the 
yield  values  for  the  shear  stresses  for  the  0°  and  ±45®  layers 
to  be  10  ksi  and  65  ksi,  respectively.  The  other  properties 
are  taken  from  table  8.  The  basic  curve  (Curve  XI)  turns  out 
to  be  a point  (fig.  52).  This  shows  that  this  notched  lam- 
inate will  fail  in  the  transverse  mode  (the  crack  running 
normal  to  the  load  direction)  prior  to  any  axial  damage. 
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This  is  based  upon  the  measured  value  (table  8)  of  the  trans- 
. T 

verse  failure  stress  a (42.4  ksi) . No  other  damage  is  pre- 
dicted to  occur  before  failure  if  this  is  the  material 
strength.  However,  this  prediction  is  based  on  assumed  shear 

yield  stress  values. 

T 

On  changing  a to  a higher  value.  Curve  I results.  It 

is  seen  that  the  0°  layers  become  inelastic  first;  the  ±45° 

layers  go  plastic  at  a higher  load  value;  and  the  laminate 

fails  i,n  the  transverse  mode  eventually.  This  behavior  can 

also  be  expected  if  the  shear  yield  stress  values  are  lower 

T 

than  the  assumed  values  and  a is  the  measured  value,  which 
may  be  the  case  for  fatigue. 

One  of  the  earlier  objectives  in  this  study  was  to  ultra- 
sonically  measure  the  degraded  properties  of  the  notched 
laminate  after  subjecting  it  to  known  fatigue  loading  condi- 
tions. This  could  not  be  achieved,  and  as  a result,  analyt- 
ical predictions  of  the  fatigue  behavior  of  notched  laminates 
(that  requires  degraded  properties  as  input  data)  were  ham- 
pered. The  fatigue  damage  growth,  though,  can  be  extracted 
from  the  ultrasonic  C-scans  (figs.  40  to  42).  It  is  seen 
that  the  extent  of  such  damage  is  small. 

The  limited  success  with  notched  laminates  restricted 
the  further  application  of  the  NDE  methodology  to  the  delam- 
ination flaw  in  composite  structures. 

DELAMINATED  BEAM  STUDY 

The  closed-form  solution  for  the  static  failure  analysis 
of  a debonded  beam  (Appendix  A)  is  plotted  in  figure  3 for 
two  locations  of  the  delamination  through  the  thickness  of 
the  midplane- symmetric  beam.  It  is  seen  that  the  critical 
load  at  which  the  debond  tends  to  propagate  increases  when 
the  debond  is  located  away  from  the  midplane.  The  same  effect 
was  observed  in  the  study  on  the  criticality  of  a delamination 
in  a laminated  cylinder  (ref.  15) . On  comparing  the  results 
for  the  beam  (fig.  3)  with  those  for  the  debonded  cylinder 
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(ref.  15),  an  interesting  contrast  is  noticed.  The  debond  in 

| * 

the  cylinder  generally  tends  to  grow  in  an  unstable  fashion  j j 

for  small  crack  lengths,  followed  by  a stable  growth  and,  ■ 

subsequently,  a region  of  debond  growth  at  a ppnstan.t,  load  j j 

(ref.  15).  This  behavior  was  observed  in  an  infinitely  long  ! ■ 

cylinder  for  two  types  of  loading:  (i)  a pressure  on  the  j 1 

_ V H 

debond  surface  and  (ii)  a constant  loading  on  the  inner  sur-  j 

face  of  the  cylinder.  In  contrast,  the  beam  analysis  of  the  ji 

present  program  (Appendix  A)  predicts  an  unstable  crack  ' { i 

growth  in  every  case  (fig.  3).  The  curve  in  figure  3,  there-  j jj 

fore,  represents  the  catastrophic  crack  growth  boundary  for  jj 

the  delaminated  beam.  j 

The  results  from  the  static  tests  on  [ (04/+452/+452/04)  ] 1 

AS/3501  beam  specimens  are  shown  in  table  10.  In  this  case, 
unstable  propagation  of  delamination,  in  the  direction  predic- 
ted by  the  finite  element  analysis,  was  observed.  The  results  ;j 

(P_„)  from  table  10  and  the  corresponding  solutions  (P_  /rv') 

cr  a i 

from  the  static  analysis  (fig.  3)  yield  a measure  of  the 
specific  adhesive  surface  fracture  energy  (table  16) . Taking 

3 i 

into  consideration  the  fact  that  a simplified  beam  analysis  was 
used,  estimated  surface  energy  values  for  two  debond  sizes  are 
in  fair  agreement.  The  quantitative  discrepancy  can  be  explained 
with  the  help  of  the  elasticity  solutions  presented  in  figure 
53.  When  the  debond  size  is  large,  the  crack  tip  nearer  the 
applied  loading  is  affected  by  the  elasticity  effect  s,hown  in 
figure  53.  Since  the  debond  tends  to  propagate  into  this  direc- 
tion, the  additional  load  needed  to  overcome  the  compressive 
stresses  at  the  crack  tip  results  in  an  increase  in  the  magni-. 
tude  of  the  critical  load.  Hence  the  larger  value  of  va  for 
the  2.54  cm  debond  in  comparison  to  that  for  the  1,27  cm  debond 
in  table  16.  "'r-  '-'-v  ■ 

Figure  11  shows  the  effect  of  fatigue  on  the  debond 
growth.  A simplified  fatigue  philosophy  is  used  to  account 
for  the  spatial  degradation  of  the  stiffnesses  due  to  cyclic 
loading,  and  typical  degradations,  from  earlier  experimental 
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observations,  are  incorporated  into  the  static/fatigue  failure 
model.  As  expected,  figure  11  shows  a decrease  in  the  criti- 
cal load  (P  ) with  an  increase  in  the  number  of  fatigue  load 
cycles  (N) . The  degradation  is  relatively  insignificant  for 
almost  a million  cycles. 

Figures  48  and  49  show  the  growth  of  the  damage  area  with 
fatigue  loading  for  two  initial  debond  sizes  in  a [(0./+45o)  ] 
laminated  beam.  The  damage  areas  in  the  figures  are  measured 
off  the  ultrasonic  "C"  scans  at  the  failure  site.  The  cor- 
respondence between  the  damaged  area  measured  off  the  "C"  scan 
and  the  damage  area  addressed  in  the  analysis  needs  further 
investigation.  If  this  relationship  can  be  established,  fig- 
ures 48  and  49  define  the  path  labeled  AED  in  figure  13. 

As  explained  in  the  section  entitled  "Analytical  Method- 
ology," figure  13  presents  the  instability  curves  for  a 
[04/±45214s  laminated  beam  with  a debond  in  its  midplane. 

One  boundary  assumes  that  there  is  no  degradation  in  the 
specific  adhesive  surface  energy  due  to  cyclic  loading.  The 
other  solid  curve  represents  the  instability  boundary  when 
the  surface  energy  degenerates  at  the  same  rate  as  the  shear 
modulus.  In  the  latter  case,  a degradation  in  the  surface 
energy  leads  to  catastrophic  failure  earlier.  This  can  be 
seen  by  following  the  actual,  stable  damage  path  AED  in  the 
figure.  Surface  energy  degradation  also  causes  a reduction 
in  the  residual  strength  and  lifetime  of  the  fatigued  struc- 
tural element.  Changes  in  EF  and  EG  in  figure  13  give  a 
measure  of  'these  reductions. 

EXPLORATORY  STUDIES 

Vibrational  monitoring  was  explored  as  a candidate  NOT  to 
measure  flaw  size  and  to  detect  flaw  growth.  Correlation  be- 
tween analytical  predictions  and  experimental  Measurements  was 
fair  (table  7).  But  further  improvement  is  necessary  before 
vibrational  monitoring  can  be  employed  as  an  NOT.  A higher 
order  analysis  may  be  required  to  achieve  this  objective. 
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Figures  32  and  33  trace  the  load-deflection  curves  for 
a t (O4/+452/+452/O4) s]  laminated  beam  with  and  without  an 
interlaminar  debond.  These  results,  while  they  do  not  form  an 
integral  part  of  the  NDE  logic,  serve  as  a potential  NDT  to 
estimate  stiffness  degradation  due  to  fatigue  loading.  While 
the  unflawed  beam  failed  in  flexure,  the  delaminated  beam 
failed  in  shear,  manifested  by  a catastrophic  growth  of  the 
debond.  The  slope  of  the  load-deflection  curve  gives  a measure 
of  the  beam  stiffness.  A comparison  of  the  slopes  of  figures 
32  and  33  shows  a fifty  percent  reduction  in  the  beam  stiffness 
for  an  a/L  of  1/3.  Similar  curves,  if  obtained  for  fatigued 
specimens,  provide  an  alternative  method  to  estimate  the 
stiffness  degradation  due  to  fatigue  loading. 
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CONCLUSIONS 

A logic  for  an  NDE  methodology  to  evaluate  the  integrity 
of  structural  elements  was  explored  under  the  premise  that 
the  goal  is  quantitative  measurement  of  residual  strength  and 
lifetime  of  flawed  components.  The  methodology  v>as  based  on: 
a static  failure  analysis  to  predict  flaw  criticality;  a sim- 
plified fatigue  philosophy  to  predict  flaw  growth;  an  experi- 
mental program  to  support  and  verify  the  analysis;  and  an  NDE 
technique  to  use  NDI  measurements  in  conjunction  with  the 
analysis  to  determine  flaw  criticality. 

A slit  notch  through  the  thickness  of  a laminate  and  an 
interlaminar  debond  in  a layered  beam  were  identified  as  crit- 
ical flaws  of  interest.  An  existent  static/fatigue  failure 
model  for  the  slit  notch  study  was  used  and  related  experi- 
ments carried  out.  Correlation  between  analysis  and  experi- 
ments was  hampered  by  many  factors.  Most  of  the  observed 
damages  were  due  to  free  edge  effects  that  are  not  accounted 
for  in  the  analysis,  and  degraded  properties  of  the  fatigued 
specimens  could  not  be  obtained  in  enough  detail. 

A static  failure  analysis  was  developed  to  predict  the 
critical  loads  for  the  growth  of  debonds  in  layered  beams 
within  the  limitations  of  a Timoshenko  beam  analysis.  A sim- 
plified fatigue  philosophy  was  adapted  to  yield  a static/ 
fatigue  failure  model  for  debonded  beams  to  establish  insta- 
bility boundaries.  The  varying  levels  of  degradation  in  the 
material  properties  with  a growth  in  the  debond  size  were 
approximated  by  those  corresponding  to  the  initial  debond 
size.  The  restrictions  involved  in  obtaining  these  bounds 
should  be  relaxed  in  future  investigations. 

An  experimental  program  was  conducted  to  support  the 
analysis  and  to  measure  the  specific  adhesive  surface  frac- 
ture energy  quantitatively.  The  discrepancy  in  the  measured 
Y values  for  two  different  debond  sizes  suggests  the  need  for 

d 

a more  refined  analysis  to  account  for  the  factors  that  cause 
it,  or  a change  in  the  experimental  method  of  obtaining  Pcr 
to  eliminate  unwanted  effects. 
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Vibrational  monitoring  of. the  specimens  to  detect  damage 
growth  and  a vibration  analysis  were  carried  out  for  correla- 
tion studies  and  to  explore  the  possibility  of  establishing 
a reliable  NDT.  Initial  success  in  the  experimental-analytical 
correlation  study  suggests  an  improvement  in  the  analysis  for 
better  agreement  in  the  results. 

Ultrasonic  "C"  scans  were  used  to  ensure  quality  control 
and  to  determine  damage  growth  rates.  The  correspondence 
between  the  damage  area  in  the  "C"  scans  and  the  damage  area 
addressed  in  the  analysis  is  not  obvious  and  calls  for  further 
study. 

The  study  carried  out  in  this  report  defines  the  potential 
of  a methodology  for  quantitatively  assessing  residual  proper- 
ties and  lifetime  of  flawed  composite  structural  elements. 

The  translation  of  this  approach  into  a quantified  criterion 
for  damage  tolerance  is  the  goal,  and  this  would  establish 
flaw  criticality  and  hence  the  need,  or  lack  of  it,  for  man- 
datory repair.  The  results  identify  areas  where  further  inves- 
tigation is  mandatory.  Limited  success  with  the  slit  notch 
study  suggests  that  effort  should  focus  on  an  improved  delam- 
ination study  to  establish  the  proposed  NDE  methodology  as  a 
workable  design  tool.  A higher  order  analysis  for  this  inter- 
laminar defect  study  and  improved  experimental  methods  for 
defining  damage  growth  are  required.  Also,  the  validity  of 
this  methodology  and  the  definition  of  a procedure  for  its 
application  to  real  structures  have  to  be  established. 
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Table  1.  Types  of  Flaws  and  Imperfections 
in  Composite  Structures 


1.  THROUGK-THE -THICKNESS  FLAWS 

a.  Slit  Notch 

b.  Openings  or  Cutouts  (circular/  elliptical/ 
rectangular /arbitrary  shape) 

c.  Bolted  Joints,  Interference-Fit.  Fasteners 

2 . SURFACE  FLAWS 

a.  Region  Of  Broken  Fibers 

b.  Breaks  Along  The  Length  Of  A Fiber 

c.  Cracks  Parallel  To  The  Fibers 

3.  SUB-SURFACE  FLAWS 

a . Interlaminar  Damage 

b.  Cracks  Parallel  To  The  Fibers 

c.  Regions  Of  Broken  Fibers 

d.  Breaks  Along  The  Length  Of  A Fiber 

e.  Matrix  Porosity 

4 . OTHER  IMPERFECTIONS 

a.  Variation  Of  Local  Fiber  Orientation  And 

Spacing 

b.  Local  Fiber  Microbuckling  (in  the  vicinity  of 

bolt  bearing  surface) 
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Table  2.  Static  Test  Data  for  t <04/i452/+452/04) gl g AS/3501  Laminates 
With  and  Without  Interlaminar  De.' aminations 
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Table  3 . 

Input  Data 

for  the  Finite  Element 

Model 

Layers 

Layers 

Property 

1 and  5 

2 and  4 

Layer 

Lay-up 

I0]16 

t±458/°i6] 

t±4516 

Thickness 

(mm) 

2.235 

4.470 

4.470 

Ex  (GPa) 

127.550 

66.360 

19.150 

Ey  - Ez  (GPa) 

10.070 

20.570 

19.150 

G (GPa) 

5.520 

17.130 

33.100 

V)  = V)  = \) 


n Tin 


0.201 


0.736 


Table  4 . Vibration  Results  for  a t s^s  AS/-*501  Laminate 

with  No  Interlaminar  Defect  (S=0.5) 
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Table  7. 


Frequency 
(Hz ) 


n.j/L  = 0.2 


L = 12.7cm 


Comparison  of  Experimental  and  Analytical 
Frequencies  of  Vibration  of  [(0./+45,)  ] 
AS/3501  Laminates  With  and  Without  Inter- 
laminar Delaminations 


a/L  m Q 

Analysis  Experiment 


a/L  - 0.2* 

Analysis  Experiment 


“l 

313 

342 

255 

340 

« "****'■•  ; 

<°2 

1958 

2049 

1653 

1961 

w3 

5483 

5407 

4635 

5351 

“4 

10745 

10343 

8812 

9603 

►47- 


. ■ * , 


....  • ~r-r- ..  ■ . , 

-vt> a- -JV:  <, . ■■■; 
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Summary  of  Static  Tension  Test  Data  for  Various  AS/3501  Laminates 
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Table  10. 


Static  Test  Results  for  [ (04/+452/+452/04) g] g 
AS/3501  Laminates  With  and  Without  Interlaminar 
Delaminations 


Width 


Thickness 


Debond  size** 


Critical 

Load 

p *** 


cr 


Specimen* 

(cm) 

(cm) 

(cm) 

(N) 

1 

2.545 

0.902 

1.27  and  2.54**** 

3139.0 

2 

2.543 

*• 

II 

3090.0 

3 

2.545 

II 

II 

2894.0 

Average: 

3041.0 

1-1 

2.535 

0.907 

2.54 

3188.5 

1-2 

2.558 

0.914 

•1 

3384.5 

1-3 

2.568 

0.887 

U 

2305.5 

1-4 

2.532 

0,881 

II 

2599.5 

Average: 

2869.5 

1/2-1 

2.545 

0.907 

1.27 

4660.0 

1/2-2 

•1 

0.914 

II 

4512.5 

1/2-3 

2.558 

II 

M 

4782.5 

Average : 

4651.5 

0-1 

2.543 

0.912 

o.o 

4856.0 (FF) 

0-2 

II 

0.884 

M 

4463. 5(FF) 

0-3 

2.548 

0.917 

II 

4610. 5 (FF) 

Average: 

4643.5 

♦Total  length  = 17.78cm?  Test  section  * 15.24cm 
**Debond  was  centered  in  the  test  section 
** ‘Catastrophic  propagation  of  delamination 
♦♦♦♦One  size  on  one  half  of  the  specimen. 

* * * * denotes  f lexural 
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•Total  length  of  specimen  « 30.48cm;  Test  Section  Length  = 22.86cm 
••Notch  sise  = 0.9525cm. 
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Table  16 . 


Computed  y&  Values  for  Two  Delamination  Sizes 


Specimen* 

a/L 

Pcr 

(lb) 

From  Analysis 

(lb. in) 1/2 

Ya 

(lb/ in) 

1-1 

1/3 

716.5 

456.4 

2.46 

1-2 

fl 

760.5 

II 

2.78 

1-3 

<1 

518.0 

II 

1.29 

1-4 

U 

584.0 

M 

1.64 

Average : 

2.04 

1/2-1 

1/6 

1047.0 

91,2.8 

1.32 

1/2-2 

U 

1014.0 

M 

1.23 

1/2  3 

M 

1075.0 

«a 

1.39 

Averages  1.31 


*Debond  is  in  the  midplane  and  is  centered  in  the  test 
section  of  length  - 7.62  cm. 
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Figure  1.  Minimechanics  Model  for  a Notched  Laminate 
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Figure  4.  Shear  Stress  Distribution  in  Orthotropic  Elastoplastic 
Short  Beam  Under  Three-Point  Bending  at  Four  Cross 
Sections:  A,  B,  C & D.  (Reproduced  from  C.  A.  Berg/ 
et  &1. / "Analysis  of  Short  Beam  Bending  of  Fiber 
Reinforced  Composites",  ASTM-STP  .4971; 
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Element  Numbering  for  the  Finite  Element  Model 
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At  the  Crack  Tip  away  from  the 
Support  (x,  » A.  + a) 

At  the  Crack  Tip  nearer 
/ — the  Support  (x  * ) 


Y I016/±458]2S  AS/3501 


Laminate 


a/L  * 0,2  (L“  11.43  cm) 

A^/L  * 0.4 

Debond  in  the  midplane  (y*0) 


Figure  10 . Centroidal  Transverse  Shear  Stress  Distributions  at 
the  Crack  Tips  for  a [0X6/±458l2s  AS/3501  Lanjinated 
Beam  from  the  Finite  Blement  Analysis 
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[(04/+452/+452/04)g]s  Laminate 

Oelamination  Between  Plies 
32  and  33 

ix/L  = 1/3 

dy/dN  = -0.21E-7 
dG12/dN  = 0.21E-6 

dE11/dN  = dE22/dN  = 0 

A corresponds  to  P_tafcie//y_ 
for  a/L=0 .2.  static  v a 

B is  the  corresponding  maximum 
fatigue  load. for  S/n=2/3. 

Nf  is  the  number  of  cycles  to 
l failure  if  a/L  remains 

\ constant  (0.2). 


• N*l  cycle 


c§c!8s 

I V 

7 1 N* 

N=10  cyclejs  


Figure  11,  Effect  of  Fatigue  Cycles  on  the  Critical  Delamination 
Load  for  a l S1  s Laminate 


NADC-76228-30 


Figure  13.  Catastrophic  Crack  Growth  Boundary  for  an 
Initial  a/L  * 0.2  and  S/n  “ 2/3  for  a 

l (°4/i452/5r'J52/04)  s1  e Lailttinate 
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Figure  17  » Typical  Stress-Strain  Curve  for  a [0]  AS/3501  Laminate 


Figure  20.  Typical  Stress-Strain  Curve  for  a [±45]  AS/3501  Laminate 
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Figure  29.  [45/0/-45/0] s AS/3501  Specimens  used  in  the 

Slit  Notch  Static  Tests 
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{ (G4/+452/+452/04)s)s  Laminate 
Specimen  1-2 

Debond  length  = 2.54cm  / 


Catastrophic  Propagation 
of  Delamination 


1 * 0 

1.27 

2.54  3.81  5.08 

6.35  j 

! 

! 

? 

K 

i 

Tip  Deflection  (mm) 

f Figure  33. 

Variation 

of  the  l (04/±452/+452/04) 

i 

* 

sls  AS/3501  ] 

Laminated 

Beam  Tip  Deflection  with 

the  Applied  j 

' 

Tip  Load 
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after  Fabrication 
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Figure  37.  Accumulation  of  Edge  Damage  after  5x10  and  10  Cycles 
in  [45/0/-45/0]  AS/3501  Specimens  Fatigued  at  30  Hertz 

with  S = 0.67  and  R = 0.1 
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Figure  40.  "C"  Scans  of  Fatigued,  Notched  [45/0/-45/0] g AS/3501 

Specimens  at  a Gate  Frequency  of  4.2  MHz 





Figure  41.  "CH  Scans  of  Fatigued,  Notched  [45/0/-45/0]  AS/3501  Specimens 

at  a Gate  Frequency  of  7.2  MHz  s 
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Figure  43.  Residual  Strength  and  Residual  Modulus  in  a 
Notched  [45/0/-45/0]  AS/3501  Laminate  when 

S = 0.67  s 
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Figure  44.  Residual  Strength  and  Residual  Modulus  in  a 
Notched  145/0/-45/0]  AS/3501  Laminate  when 

S = 0.8  s 
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Figure  48.  Growth  of  Damage  Area  at  the  Failure  Site  with  N for  a [(04/+45 
Laminate  with  a 1.27cm  Debond 
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Figure  49.  Growth  of  Damage  Area  at  the  Failure  Site  with  N for 
Laminate  with  a 2.54cm  Debond 
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■ Initiation 
of  inplane 
inelasticity 
in  the  0° 
layers 


Initiation 
of  inplane 
inelasticity 
in  the  *45° 
layers 


Transverse 
failure  of 
the  notched 
laminate 


a - 100  ksi  for  case  I 
T 

a =42.4  ksi  for  case  II  (actual  value) 


Transverse  failure 
of  the  notched 
laminate 


[45/0/-45/0]  Laminate 
s 


Total  Damage  Length  (mm) 


Figure  52 . A Static  Failure  Analysis  of  the  Notched 

[45/0/-45/0]  Laminate  for  Two  Values  of  the 
Transverse  “Failure  Stress,  aT 
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Figure  53.  Lines  of  Constant  Stresses  for  Different 

Elastic,  Orthotropic  Half  Planes  (Reproduced 
from  "Anisotropic  Plates"  by  S.  G.  Lekhnitskii) 


-109- 


NADC-76228-30 


APPENDIX  A 

CRITICALITY  OF  INTERLAMINAR  DELAMINATIONS 
IN  LAMINATED  BEAMS 

A simplified  analysis  is  presented  here  to  estimate  the 
criticality  of  an  interlaminar  delamination  in  a laminated 
be*n>.  The  debond  may  be  a birth  defect  that  occurs  during  the 
fabrication  process  or  a service  defect  that  is  a result  of 
cyclic  loading,  impact,  etc.  Considering  a typical  practical 
application  to  be  an  aircraft  wing,  the  structural  element  is 
modeled  as  a cantilevered  laminated  beam.  The  delamination, 
through  the  width,  is  assumed  to  exist  at  any  axial  location 
along  the  length  of  the  beam,  and  between  any  two  laminae 
through  the  thickness  (fig.  2 ) . The  flawed  beam  is  then 
analyzed  as  four  different  beams  that  are  joined  together, 
with  appropriate  boundary  and  matching  conditions  imposed  to 
realize  the  true  response  of  the  debonded  beam. 

The  two  Euler  equations  for  a beam,  using  Timoshenko 
beam  theory  to  include  shear  deformation  effects,  are  (ref. 19) 


~ (El  ^-) 
3 x 1 3x; 


+ * - *>  - 3t 


i 


= o 


(Al) 


and 


m3 ' 


3 1 


i-h 


$)  3 - p = 0 


(A2) 


The  deflected  beam  geometry,  its  coordinates,  the  loading,  and 
transverse  shear  force  conventions  are  shown  in  figure  12 . 
ij>  is  the  slope  of  the  deflection  curve  neglecting  shear  de- 
formation and  8 is  the  angle  of  shear  at  the  neutral  axis. 

The  total  slope  is 

U - ip  + 3 <A3) 

E is  the  Young's  modulus  of  the  beam  material,  and  I is  the 
cross-sectional  moment  of  inertia,  m is  the  mass  per  unit 
length  of  the  beam  and  Ip  is  the  mass  moment  of  inertia 
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about  the  neutral  axis  per  unit  length  of  the  beam.  p(x,t)  is 
the  instantaneous  lateral  load  intensity  per  unit  length  at 
time  t , and 

k = k'  AG  (A4) 

where  A is  the  cross-sectional  area  of  the  beam,  G is  the  shear 
modulus  of  elasticity,  and  k'  is  a numerical  constant  depending 
on  the  geometry  of  the  cross  section.  For  a rectangular  cross- 
section,  k'  = 2/3.  For  a uniform  beam,  m = pA  and 

Ip  * pAR2,  (A5) 

where  p is  the  mass  density  of  the  beam  material  and  R is  the 
radius  of  gyration  of  the  cross  section  about  an  axis  perpen- 
dicular to  the  plane  of  motion  and  through  the  neutral  axis. 

The  static  analysis  considers  a tip-loaded  uniform  canti- 
lever beam  (fig.  2 ) . The  Euler  equations,  for  this  case, 
reduce  to: 

2 

El  2-1  + k (§£  - ip)  =0  (A6) 

dx  dx 

& [k<f£>  - *>1  - o (A7) 

The  bending  moment,  M,  and  the  transverse  shear  force,  Q,  at 
any  section  of  the  beam,  are  given  by: 

M = EI  fx'  (A8) 

and 

Q - -k  - ip)  (A9) 

The  deflections,  slopes  and  material  properties  of  the 
four  different  beams  are  identified  by  the  use  of  proper  sub- 
scripts. The  axial  coordinate  x^  for  the  ith  beam  has  its 
origin  at  the  left  end  of  that  beam  (fig.  2) . Equations 
(A6)  and  (A7)  may  be  combined  to  give  a third  order,  ordinary 
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ordinary  differential  equation  for  y^(x).  This  results  in 
twelve  constants  of  integration  when  all  the  four  beam  elements 
are  considered.  In  addition,  the  transverse  shear  force  divi- 
sion between  the  elements  above  and  below  the  crack  of  length 
a (fig.  2 ) leads  to  an  extra  unknown  to  be  solved  for.  If  the 
applied  tip  load  is  P and  the  beam  element  above  the  crack 
(i=3)  carries  a transverse  shear  force  of  £,  the  element  below 
the  crack  (1=4)  carries  a shear  force  of  (P-i*)  . 

The  thirteen  boundary  and  matching  conditions  used  in  ob- 
taining the  deflection  solution  at  any  axial  location  of  the 
debonded  beam  are: 


yl  *xl=03==0,  or,  yx=0  at  x1=0 
iK(x.,=0)  = 0 
M2 (x2=£2)  = 0 

M1(X1~£1)  = ~P<*2+a) 

W0)  = ^(x^) 

i^3(x3=0)  = ^ (x^=0) 

y3(x3=0)  = y^x^) 
y3(x3=0)  = y4(x4=0) 
y2 (x2=0)  = y3 (x3=a) 
y2(x2=0)  = Y4(x4=a) 

\p2  (x2=0)  = ipj  (x3=a) 

^2(x2=0)  = ^4 (x4=a) 

Ml<xrV  = M3(x3=03  + Vxr0)  + Maxial 


(A10) 


(All) 


(A12) 


(A13) 


(A14) 


(A15) 


(A16) 


(A17) 


(A18) 


(A19) 


(A20) 


(A21) 


(A22) 


Maxial  lending  moment  due  to  the  couple  created  by 

the  horizontal  or  axial  force  reactions  H3  and  H4  in  beam 
elements  3 and  4,  respectively.  The  axial  forces  are  a result 
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of  an  elongation  of  element  3 and  a corresponding  contraction 
of  element  4 in  order  to  ensure  compatibility  of  displacements 
and  slopes.  The  axial  elongation  in  element  3,  u3,  is  given  by: 

d4 

u3  = -Aip  + uq  (A23) 

where 


Ai|j  = ^2  (x2=0)  “ (A24) 

d4/2  - d/2  - d3/2  (A25) 


and  uQ  is  the  axial  rigid  body  displacement  in  elements  3 and  4. 
d = dj^  = d2  is  the  total  thickness  of  the  debonded  beam,  and  d3 
and  d^  are  the  thicknesses  of  elements  3 and  4,  respectively. 

It  follows  that  the  axial  contraction,  u4,  in  element  4 is: 

d3 

u4  = Uq  + Alp  — (A26 ) 


Since 


Ha 

u3  = on 


H4a 

and  u4  = aTe" 


(A27) 


3 3 4 4 

the  axial  force  equilibrium  equation  at  the  crack  tip  locations, 


h3  + h4  = 0, 


(A28) 


yields  a solution  for  the  rigid  body  displacement: 


Aip 

2(a3e3;a4e4) 


(A3E3d4-A4E4d3) 


(A29) 


if  b is  the  width  of  the  beam,  A3=bd3  and  A4=bd4>  This  simpli- 
fies equation  (A29)  to: 


u 


o 


d3d4(E3-E4) 
2(d3E3+d4V  » 


(A30) 


An  interesting  inference  from  the  above  equation  is  that 
a rigid  body  translation  takes  place  only  when  E3  is  not  equal 
to  E4.  If  the  beam  material  is  isotropic  or  if  the  delamination 
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in  a midplane-symmetric,  laminated,  anisotropic  beam  is  in  the 
midplane,  there  will  be  no  rigid  body  axial  displacements. 

When  the  debond  occurs  between  two  layers  such  that  E3  and  E 4 
have  different  magnitudes,  a rigid  body  translation  in  elements 
3 and  4 is  essential  for  displacement  and  slope  compatibility. 

With  the  above  expressions  for  the  axial  forces  in  ele- 
ments 3 and  4,  half  the  total  beam  thickness  apart,  Maxial  roay 
be  written  as: 


i i 


“axial  * ' H3d/2  “ -A3V3<l/2a 


(All) 


Using  the  boundary  and  matching  conditions,  (A1Q)  through 
(A22),  in  the  general  solutions  to  the  beam  equations  (A6)  and 
(A7)  for  all  the  four  elements,  the  deflection  function  for  the 
tip-loaded  delaminated  beam  is  obtained.  Computing  the  tip 
deflection  from  this  solution,  Clapeyron's  theorem  is  used  to 
obtain  the  total  strain  energy,  U,  in  the  beam: 


— 2 ^ ^2  ^x2-^2^ * 


(A32) 


2U/P  = - 


y2(x2=V 


(A33) 


Carrying  out  the  aforementioned  substitutions,  one  obtains: 


2 

3TeT3 


C 3 Hr) 

Hr. + r— 

2 P k2 


(A34) 


where 


C2  ^0  a (^2*^"^  *^) 


(A35) 


- f a3 

p l6(EI) 


_ _ a_i  S S 

3 kj]  P + ~^+~d  ~ 


(A36) 


/yi  aaii>sp\3>a!g>L  !•  »gawB  y - - swmswbswwp  tp9 
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Substituting  equation  (A34)  into  equation  (A45) : 


2 (3U/da) _ 4bYa  _ i2  dl2  C2  dZ2  „ d(C2/p) 

^2  ~P  (El)  2 da  " P da  “ Z2  da 

cr  cr 


d<C3/P>  _1  ^2 

da  + k2  da 


pcr^s  Ya  = 2Jb^: 


(A46) 


(A47) 


Pcr  is  the  magnitude  of  the  critical  load  at  which  propagation 
of  delamination  is  imminent.  The  derivatives  that  appear  in 
the  expression  for  <j>3,  equation  (A46),are: 


d(C2/P) 


d(C8/P) 


d(C3/P) 


a £ a2  d$/P)  C7  d«VP) 

(EI)3  P 2 (El)  3 da  P a da 


_1  £ a d (%/P)  a2  | a3 

k3  P k3  da  2 (El)  3 P 6(EI)3 

c7  a2  <*(C,/P)  C„  d(C„/P)  d (C_> 


*aa±  + + 2i™7/p)  , c8  , ad(Vp)  ,d<vp> 

A»  P 9 A^k  ▼ IT  + a T~_ + 


(A48) 


(A49) 
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d(C?/P) 

da 


[dJl2/da+0.5] 


2 (El ) 


(A50) 


d(C8/P) 

da 


(L-£1)  dJ^ 


(El)  ^ da 


(A51) 


d(Cg/P) 

da 


(L-0.5Jl1) 

(El), 


^1 

'1  da 


JL 

k1  da 


(A52) 


d(^/P)  . _ _ _ a 

da  6 (El)  4(3^2)  11  P U+(EI)  3^J 


(A53) 


It  is  seen  from  equation  (A47)  that  the  critical  value 
of  the  tip  load  is  a function  of  the  location  and  length  of 
delamination,  the  material  and  geometric  properties  of  the 
beam,  and  the  specific  adhesive  surface  fracture  energy  for 
the  given  adherend-adhesivc  combination . The  propagation  of 
the  delamination,  along  the  interface,  may  take  place  in  any 
of  the  following  ways  depending  on  the  material  and  geometric 
data : 

(i)  propagation  of  the  delamination  into  element  1; 
in  this  case,  dS-^-da  and  di2=0; 

(ii)  propagation  of  the  delamination  into  element  2; 
in  this  case,  d£2=-da  and  dJl^O; 

(iii)  propagation  of  the  delamination  into  elements  1 
and  2 simultaneously;  in  this  case,  di^di^-- 

The  critical  load  for  each  case  is  obtained  by  substituting 
the  appropriate  values  of  di^/da  and  d&2/da  into  equation 
(A47) . The  direction  in  which  the  delamination  propagates  is 
governed  by  the  least  of  the  three  Pcr  values. 

The  variation  of  with  the  crack  length  a is 

plotted  for  different  sets  of  beam  data.  One  experimental 
data  point  on  this  curve  would  suffice  to  evaluate  Ya* 
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APPENDIX  B 

VIBRATION  ANALYSIS  OF  LAYERED  BEAMS 
WITH  DELAMINATIONS 

The  vibration  analysis  of  a laminated  beam  is  carried  out 
by  eliminating  ip  from  equations  (Al)  and  (A2)  to  obtain  the 
well-known  Timoshenko  equation  for  lateral  vibration  of  pris- 
matic beams  (ref.  19  ) : 


„„  3%  . _ 3 Y ,,  , Elm,  3 y , m 3 y _ 

El  — \ + m — \ - (I  + -j— ) — y-S-  + ID  v — 4 P + 

3x4  3t  p k 3x  3t  p k 3t 

12  2 
L_£  - El  r_p 

k a?  k 3x 


(Bl) 


Setting  p=0  in  the  above  equation,  the  Timoshenko  equa- 
tion for  lateral,  free  vibration  is  obtained: 


34V  . _ 32v  ,,  , Elm,  34V  . » m 34y  _ „ 

EI  — T + m — J “ + "XT-)  — + Jo  k — t " 0 

3x  3t  p k 3x  3t  p k 3t 


Assuming 


iwt 


y (x,t)  = y (x)  e 
equation  (B2)  is  rewritten  as: 
,4 


EI  - m</y  + (I  + liS)  + I | u*y  = o 


m 4. 


dx 


dx* 


(B2) 


(B3) 


(B4) 


or. 


+ W2  (— y + — y)  - -y-yi  y « 0 

dx4  C2  C2  dx2  C2c2  C2R 

0 4 0 4 0 


(B5) 


where 


C2  *»  E/p,  C2  ® k*G/p  and  R2  * I/A 
A general  form  of  the  solution  to  equation  (B5)  is: 


(B6) 


y(x)  » B sin  ax  + C cos  ax  + D sinh  &x  + E cosh  8x  (B7) 

B-l 
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where  B,C,D  and  E are  constants.  If 


2 , i . 1. 

“ + nr> 

o Cq 


(B8) 


and 


u,4 

TT 

o q 


co2 

2T? 


CoR 


(B9) 


a and  3 are  the  solutions  to  the  following  quartic  equations: 

(BIO) 
(Bll) 


4 2 , . 

a - ya  +6=0 


6 4 + y82  +6=0 


The  four  constants  in  equation  (B7)  are  evaluated  through 
the  application  of  the  proper  boundary  conditions.  As  dis- 
cussed in  Appendix  A,  the  delaminated  beam  is  treated  analyt- 
ically as  four  beam  elements  joined  together  with  the  appro- 
priate boundary  and  matching  conditions  (fig.  2 ) . Sub- 
scripts (i=l,4)  on  the  various  quantities  identify  the 
elements  they  correspond  to.  Sixteen  conditions  are  specified 
to  evaluate  the  free  vibration  f requenci  'ss  that  yield  non- 
trivial solutions  for  the  constants  B^,  and  E^: 


yi<*i*o)  * o 

(B12) 

^l (xj*Q)  = o 

(B13) 

Q2(X2*V  = 0 

(B14) 

M2(x2*V  ‘ 0 

(B15) 

y1(x1»£1)  = y3(x3»0) 

(B16) 

yl(xl=“£l)  * y4  (x4*0) 

(BIT) 

ij/1(x1=a1)  = tjj3(x3«0) 

(B18) 
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^l(xl=il)  ■ ^4(x4*0) 

(B19) 

y3(x3=a)  = y2(x2=0) 

(B20) 

y4(x4=a)  * y2(x2=0) 

(B21) 

i]>3(x3=a)  * \p2  (x2=0) 

(B22) 

*4(x4=a)  = «2(x2=0) 

(B23) 

Q1(x1-Jl1)  = Q3(x3=0)  + Q4(x4=0) 

(B24) 

VW  “ M3(x3=0)  + M4<x4=0)  + Maxial 

C2/'X2*0)  = Q3(x3=a)  + Q4(x4=a5 
M2(x2=0)  » M3<x3~a)  + M4(x4=a)  + Maxial 


(B25) 

(B26) 

(B27) 


In  imposing  boundary  conditions  that  involve  the  bending 
slopes,  the  general  form  of  the  bending  slope  solution  (ref. 20) 
is  used : 

2 2 

!j)(x)  = -k-,(^E-a£  * g ID  cosh  Bx  + E sinh  Bx]  + 


B 

(k ' G/E) p2  - a2  , _ . _ . 

— ^ l-B  cos  ax  + C sxn  ax] 


(B28) 


where 


2 ..2  y_2 
P = w /CQ 


(329 ) 


Equations  (Bl2)  to  (B27)  and  the  solutions  (B7)  and  (B28)  yield 
sixteen  homogeneous  equations  for  the  four  sets  of  constants 
Bi*  Ci'  Di  Ei  • These  equations,  when  written  in 

matrix  form,  yield  non-trivial  solutions  for  the  constants  when 
the  determinant  of  the  matrix  vanishes . 

The  free  vibration  frequencies  are  obtained  numerically  by' 
computing  the  determinant  of  the  16x16  matrix  for  different 
values  of  frequencies  and  finding  the  frequencies  corresponding 
to  a null  determinant  by  interpolation.  These  frequencies  may 
be  normalized  with  respect  to  the  fundamental  frequency  of 
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